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ABSTRACT 

The  kinetics  of  the  reversible  binding  of  cyanide  by 
ferric  horseradish  peroxidase  was  studied  at  25*0°  and  an 
ionic  strength  equal  to  0°11  over  the  pH  range  4®  2  to  11 ® 3 
by  means  of  a  stopped- f low  apparatus «  Analysis  of  the  pH 
dependence  of  the  bimolecular  rate  constants  is  consistent 
with  the  presence  of  three  heme-linked  ionizable  groups  on 
peroxidase  with  pK  values  of  4ol9  6o4  and  10® 8 o  It  is  not 
possible  to  distinguish  whether  HCN  or  CN"  is  the  attacking 
species  but  the  bound  form  of  the  ligand  appears  to  be  the 
cyanide  ion. 

The  optical  rotatory  dispersion  of  horseradish  peroxi- 
dase  and  its  cyanide *  fluoride  and  hydroxide  complexes  was 
studied  in  the  spectral  region  from  215  mp  to  450  mp»  and 
that  of  the  azide  complex  from  350  mp  to  450  mp.  The  effect 
that  splitting  the  heme  from  the  protein  of  peroxidase  has 
on  the  optical  rotatory  dispersion  in  the  215  mp  to  450  mp 
region  was  also  studied.  Results  of  measurements  of  the 
reduced  mean  residue  rotation  at  233  mp  lead  to  the  conclu¬ 
sion  that  there  are  no  significant  changes  in  gross  protein 
conformation  upon  the  binding  of  ligands  to  peroxidase $  but 
that  the  splitting  of  the  heme  causes  a  reduction  of  the 
helical  content  of  the  protein.  Pure  peroxidase  was  esti¬ 
mated  to  have  43  °/o  ^“helical  content ,  which  was  reduced 
to  33  °/0  helical  content  when  the  heme  was  split  from 
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the  protein o  Results  of  studies  in  the  Soret  region  indi¬ 
cate  that  the  binding  of  various  ligands  does  not  cause  an 
alteration  of  the  geometry  of  the  heme  with  respect  to  the 
protein  moiety. 

The  equilibrium  constants  for  the  formation  of  the 
hydroxide  complex  of  ferric  horseradish  peroxidase  were  de¬ 
termined  spectrophotometrically  at  25°0°  for  seven  values 
of  ionic  strength  over  the  range  0.02  to  0o20  and  also  at 
I8o0°  and  35o0°  with  p.  =  0.11.  The  results ,  which  were  ana¬ 
lyzed  on  the  assumption  that  the  hydroxide  complex  is  formed 
by  the  ionization  of  a  water  molecule  in  the  sixth  coordin¬ 
ation  position  of  the  heme  ferric  iron,  lead  to  the  conclu- 
sion  that  two  isozymes  of  peroxidase  are  present  in  the 
preparations  of  horseradish  peroxidase  which  were  used.  The 
isozymes  differ  in  that  they  have  different  values  for  the 
pK  of  peroxidase-hydroxide  complex  formation.  The  efiec- 
tive  charge  on  the  aquo  species  of  the  predominant  isozyme 
is  indicated  to  be  between  +2  and  +3 ;  its  dissociation  con¬ 
stant  is  (2olT0o2)  xlO"11M  and  the  thermodynamic  parame¬ 
ters  for  its  ionization  are  AG°  =  l^o 6^0.1  kcal/mole;  — 

7o8±Oo 1  kcal/mole?  As0  =  -22o6±lo0  entropy  units  at  25°  and 
zero  ionic  strength.  Results  of  temperature- jump  relaxa¬ 
tion  kinetic  experiments  indicate  that  the  rate  constant  for 
the  combination  of  a  proton  with  the  peroxidase-hydroxiae 
complex  is  greater  than  1.3  x  lO10  M”lsec  1  at  31 °° 

The  NMR  halide  ion  probe  technique  was  shown  to  be  a 
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very  sensitive  method  for  detecting  the  binding  of  Cl  and 
Hgs+  to  proteins .  There  appears  to  be  only  a  small  amount 
of  chloride  ion  binding  to  horseradish  peroxidase  and  to 
its  apoprotein  in  basic  aqueous  solutions,  but  the  binding 
increases  markedly  as  the  pH  is  lowered,  indicating  that 
the  chloride  ion  is  binding  non-specifically  to  positively 
charged  regions  on  the  protein  <,  HgST  was  found  to  bind  to 
myoglobin  in  increasing  amounts  as  the  pH  was  raised  above 
7,  indicating  that  it  may  be  binding  to  regions  of  negative 
charge  on  the  protein o  Use  of  the  halide  ion  probe  tech“ 
nique  as  a  means  of  detection  for  HgCla  titrations  resulted 
in  the  finding  that  there  is  one  reactive  sulfhydryl  group 
in  peroxidase o  The  nature  and  number  of  the  sulfhydryl 
groups  of  hemoglobin  and  catalase  wer e  also  investigated o 
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CHAPTER  1.  INTRODUCTION 


The  organization  of  this  thesis  reflects  the  fact  that 
a  number  of  different  projects  were  undertaken  in  the  course 
of  the  author's  research.  Hence,  the  thesis  consists  of 
several  chapters,  each  of  which  is  meant  to  be  complete  in 
itself  with  regard  to  an  introduction,  experimental  details, 
and  discussion  of  the  results  for  each  research  problem. 

However,  since  the  major  part  of  the  research  was  car¬ 
ried  out  on  the  enzyme  horseradish  peroxidase  (abbreviated 
to  HRP),  the  following  chapter  is  included  to  acquaint  the 
reader  with  the  nature  and  physical  properties  of  HRP. 

General  Introduction  to  the  Nature  of  Peroxidase.  Its 

Occurrence  and  Isolation 

Peroxidase  is  an  enzyme  which  catalyzes  oxidations  in 
the  presence  of  hydrogen  peroxide  and  has  been  given  the 
following  designation  by  the  Enzyme  Commission  of  the  Inter¬ 
national  Union  of  Biochemistry1:  1.11. 1.7  Donor:  H202 
Oxidoreductase.  It  is  much  less  specific  in  its  action  than 
many  enzymes,  with  known  substrates  including  leuco-dyes, 
phenols,  indole,  aromatic  amines,  ascorbic  acid,  and  some 
inorganic  ions  such  as  iodide.  A  discussion  of  specific 
oxidative  reactions  involving  peroxidase  and  a  review  of 
the  theories  of  the  mechanism  of  its  enzyme  action  are 
found  elsewhere3,  and  will  not  be  discussed  further  here. 


Peroxidases  are  present  in  human  saliva,  leucocytes, 
animal  thyroid  tissue,  milk,  and  many  higher  plants.  Hence¬ 
forth  the  material  in  this  chapter  will  specifically  refer 
to  the  peroxidase- obtained  from  its  most  common  source, 
horseradish  roots. 

The  first  isolation  of  HRP ,  in  partially  purified  form, 
was  carried  out  by  Bach  and  Chodat  in  19033*  It  was 
crystallized  by  Theorell  in  1942  using  electrophoretic  tech¬ 
niques4.  Other,  simpler  methods  of  preparing  pure  HRP  have 
been  described  by  Keilin  and  Hartree6  and  by  Paul6.  HRP  is 
now  available  commercially  in  various  degrees  of  purity 
from  a  number  of  firms. 

HRP  is  known  to  consist  of  a  number  of  components  which 

vary  in  peroxidatic  activity,  amino  acid  composition,  and 

spectral  properties.  The  first  evidence  of  these  multiple 

components,  or  isozymes,  was  found  by  Theorell,  who  isolated 

a  form  of  HRP  which  he  called  peroxidase  I  or  paraperoxidase 
* 

Since  that  time,  studies  on  peroxidase  isozymes  have  been 
made  by  Jermyn  and  Thomas8  and  by  Klapper  and  Hackett®. 
Recently,  a  thorough  and  detailed  study  of  the  physical  and 
chemical  properties  of  the  isozymes  was  reported  by  Shannon, 
Kay  and  Lew10'  11 . 

Structure  and  Physical  Properties 

HRP  has  a  molecular  weight  of  40,20012,  and  an  iso¬ 
electric  point  of  7.213.  It  is  thermally  stable  for  15  min. 


. 


-3- 


at  63 0 »  and  for  weeks  at  room  temperature14.  Elemental 
analysis  shows  HRP  to  contain  47*0  °/o  carbon,  7*25  °/o 
hydrogen,  13»2  °/o  nitrogen,  0.43  °/o  sulfur,  0.127  °/o 
iron,  and  32  °/o  oxygen13.  It  has  been  found  that  about  18 
weight  percent  of  HRP  is  carbohydrate  material10’13,  and 
although  the  nature  of  the  linkage  to  the  protein  is  not 
known,  the  composition  of  the  carbohydrate  has  been  stud¬ 
ied9'  10.  The  protein  moiety  contains  287  amino  acid  resi¬ 
dues13;  the  sequence  of  these  has  not  been  determined,  but 
qualitative  and  quantitative  amino  acid  analyses  have  been 
published  by  a  number  of  workers9 > 10> 1 5> 16 . 

HRP  is  a  member  of  the  family  of  heme  proteins,  which 
also  includes  hemoglobin,  myoglobin,  and  catalase.  The 
common  feature  of  these  proteins  is  that  they  all  contain 
as  a  prosthetic  group  protohematin  IX,  shown  in  Figure  1. 
Hemoglobin  and  catalase  contain  four  heme  groups  and  HRP 
and  myoglobin  each  contain  one  heme  group.  The  hematin 
nature  of  HRP  was  established  by  Keilin  and  Mann17,  and  the 
importance  of  the  heme  group  was  demonstrated  by  Theorell 
who  found  that  neither  the  protein  nor  the  prosthetic  group 
independently  has  peroxidatic  activity18.  In  its  natural 
form,  the  iron  in  the  heme  group  of  HRP  is  in  the  ferric 
state.  It  can  be  reduced  with  sodium  hydrosulfite,  but  is 
re-oxidized  in  the  presence  of  oxygen.  The  reduction  pot¬ 
entials  of  HRP,  determined  at  30°  by  Harbury,  were  found  to 
range  from  -0.1002  at  pH  3 <*95  to  -0.3372  at  pH  11.2519* 
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Figure  1*  Protohematin  IX 
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The  heme  group  is  bound  to  the  protein  by  several 
types  of  bonds,  with  the  nature  of  those  bonds  being  best 
known  for  myoglobin,  since  its  crystal  structure  has  been 
determined  in  detail  by  X-ray  studies.  The  heme  of  myo¬ 
globin  is  located  in  a  cleft  of  the  protein  with  the  vinyl 
groups  of  the  porphyrin  ring  buried  in  a  hydrophobic  region 
and  the  propionic  acid  groups  bonded  to  the  surface  by 
either  hydrogen  bonds  or  salt  linkages20.  It  has  been  es¬ 
tablished  that  a  histidine  residue  and  water  occupy  the 
fifth  and  sixth  coordination  positions  of  the  heme  iron, 
with  another  histidine  residue  hydrogen  bonded  to  the  water31. 
X-ray  studies  on  a  number  of  high- spin  ferric  porphyrins 
have  shown  that  in  all  cases,  the  iron  atom  is  displaced 
greater  than  0.3  A  out  of  the  plane  of  the  nitrogen  atoms, 
presumably  because  the  high- spin  ferric  iron  has  an  electron¬ 
ic  structure  which  gives  it  too  large  an  effective  size  to 
be  accommodated  in  the  plane  of  the  nitrogen  atoms33.  For 
myoglobin,  this  displacement  is  toward  the  bound  histidine30. 

To  date  there  has  been  no  X-ray  crystallographic  work 
reported  for  HRP,  so  evidence  regarding  the  nature  of  the 
heme  -  protein  forces  is  mostly  based  on  a  comparison  of 
the  properties  of  HRP  with  those  of  hemoglobin  and  myoglobin. 
There  has  been  a  considerable  amount  of  speculation  in  regard 
to  the  identity  of  the  ligands  occupying  the  fifth  and 
sixth  positions  on  the  heme  iron  of  HRP.  On  the  basis  of 
a  comparative  study  of  the  reduction  potentials,  magnetic 
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moments,  and  spectra  of  the  heme  proteins,  Brill  and 
Williams  concluded  that  HRP  is  "an  amino- carboxylate,  or 
less  probably,  an  amino-water  complex”23.  However,  Brill 
later  discussed  the  formation  of  HRP-ligand  complexes  in 
terms  of  the  ligand  replacing  water  in  the  sixth  position24, 
a  scheme  also  proposed  by  Nicholls25.  Additional  support 
for  the  belief  that  water  is  in  the  sixth  position  is  ob¬ 
tained  from  a  study  of  the  temperature  dependence  of  the 
spectra  of  HRP  at  various  values  of  pH6.  Nicholls  and 
Brill  both  disagreed  with  the  views  of  George  and  Lyster, 
who  postulated  that  the  heme  of  HRP  is  in  a  protein  crevice 
with  protein  residues  bound  in  both  the  fifth  and  sixth 
positions26.  Recently,  Brill  and  Sandberg  concluded  from  a 
difference  spectroscopy  study  in  the  210-280  njjti  region  that 
histidine  must  be  in  the  fifth  position  of  HRP27. 

Early  titration  data  from  the  work  of  Theorell  and  Paul 
implicated  the  propionic  acid  groups  of  the  porphyrin  in  the 
binding  of  the  heme  to  the  protein  of  HRP28.  This  has  been 
substantiated  by  experiments  performed  by  Maehly,  who 
showed  that  protohematin  IX-dimethylester  and  protoporphyrin 
IX-dimethylester  (no  iron  present  in  the  ring)  do  not  bind 
to  the  apoprotein  but  that  free  protoporphyrin  IX  does  bind 
to  apo-HRP29. 

A  detailed  study  of  the  splitting  of  HRP  into  heme  and 
protein  has  been  reported  in  a  series  of  papers  by  Maehly30"32 
HRP  is  stable  between  pH  4.2  and  11.4.  The  rate  of  the 
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splitting  reaction  depends  on  [H+]a  or  [0H~]2,  thus  it  in¬ 
creases  markedly  when  the  pH  is  lowered  or  raised  outside 
the  range  of  stability.  The  half-time  of  the  splitting 
reaction  is  about  one  minute  at  either  pH  2.5  or  13 •  At  low 
pH,  the  splitting  is  reversible  and  takes  place  in  four 
spectroscopically  distinguishable  steps;  at  high  pH,  only 
one  step  is  observed  and  the  splitting  is  irreversible.  The 
splitting  reaction  at  low  pH  occurs  faster  at  any  given  pH 
in  the  presence  of  ligands  which  bind  to  the  sixth  position 
of  the  iron  of  HRP.  As  the  heme  is  split  from  the  protein, 
there  is  an  accompanying  loss  of  peroxidatic  activity, 
which  is  regained  if  the  solution  is  neutralized  to  allow 
the  heme  and  protein  to  recombine. 

Ligand  Binding  to  HRP 

The  heme  proteins  perform  their  biological  functions 
by  interaction  with  ligands,  such  as  in  the  transport  of 
oxygen  by  hemoglobin.  Although  magnetic  and  spectral  evi¬ 
dence  indicates  that  the  iron  of  the  heme  is  involved,  the 
most  positive  proof  is  obtained  from  X-ray  crystallographic 
studies,  which  show  that  for  myoglobin,  the  ligand  is  found 
in  the  sixth  coordination  position  of  the  iron,  having  re¬ 
placed  the  water  molecule  normally  present  in  that  posi¬ 
tion31’  33  .  It  is  generally  accepted  that  ligands  bind  to 
the  sixth  position  for  all  other  heme  proteins  as  well. 

A  number  of  ligands  bind  to  ferri-HRP  with  accompanying 
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changes  in  the  spectral  and  magnetic  properties  of  the 
enzyme.  Figures  2  and  3>  taken  from  the  work  of  Keilin  and 
Hartree6,  show  the  absorption  spectra  of  ferri-  and  ferro- 
HRP  and  the  hydroxide,  cyanide,  azide,  and  fluoride  complex¬ 
es  of  ferri-HRP.  The  positions  of  the  absorption  peak  max¬ 
ima,  along  with  their  corresponding  molar  absorptivities 6 
are  summarized  in  Table  I,  which  also  lists  the  magnetic 
moments  of  the  various  species  as  determined  by  Theorell34. 

Table  I 

Spectral  and  Magnetic  Properties  of  HRP 


species 

and  Ferri-HRP  Ligand  Complexes 

wavelengths  (rryu)  of  absorption  maxima 
and  corresponding-  values  of  £  x  10" 4 

magnetic 

moment 

ferri-HRP 

642(0.284) 

497(1.003) 

403(9.10) 

5.45 

B.M 

f erro-HRP 

556(1.122) 

- 

437(7.73) 

5.19 

B.M 

HRP-OH 

574(0.686") 

545(0.864) 

416(8.85) 

2.66 

B.M 

HRP-F 

612(0.693) 

488(0.755) 

404(13.00) 

0 

CT\ 

0 

B.M 

HRP-CN 

538(1.020) 

- 

423(9.39) 

2.65 

B.M 

hrp-n3 

635(0.170) 

534(0.820) 

4l6 (11.40 ) 

-> 

The  ferric  iron  in  both  HRP  and  HRP-F  appears  to  be  in  a 
high-spin  form,  whereas  for  HRP-OH  and  HRP-CN  it  is  low-spin. 
In  addition  to  the  complexes  listed  in  Table  I,  ferri-HRP 
has  been  shown  to  react  with  hydrogen  sulfide,  with  the  re¬ 
sulting  complex  having  a  magnetic  moment  of  2.40  Bohr  Mag¬ 
netons34.  Nitric  oxide  forms  complexes  with  both  ferri-  and 
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gure  3.  Absolute  absorption  spectra  (molar  absorptivity  vs.  wave 
■ngth)  of  f erri-HRP ,  its  hydroxide  complex  and  of  ferro-HRP. 
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ferro-HRP;  stable  complexes  are  also  formed  by  ferro-HRP 
with  carbon  monoxide  35  and  with  cyanide36. 

Early  quantitative  studies  on  ligand  binding  by  ferri- 
HRP  included  the  investigations  by  Theorell34  and  by 
Chance37  on  the  effect  of  pH  on  the  equilibrium  constant  of 
the  reaction  between  HRP  and  fluoride,  and  Chance’s  kinetic 
and  equilibrium  study  of  the  reaction  of  HRP  with  cyanide 
at  two  values  of  pH38.  There  has  been  a  great  deal  of  spec¬ 
ulation  as  to  whether  the  attacking  and/or  bound  form  of  the 
ligand  is  the  neutral  species  or  the  aniona4"a6>  39;  4°.  In 
1967»  Dunford  and  Alberty  reported  a  temperature- jump  kinetic 
study  of  fluoride  binding  to  HRP  and  showed  that  the  var¬ 
iation  of  the  apparent  bimolecular  rate  constants  with  pH 
could  be  explained  by  a  mechanism  involving  the  fluoride  ion 
attacking  two  forms  of  HRP41 •  The  two  forms  differ  only  by 
the  presence  of  a  proton  on  a  heme-linked  ionizable  group 
with  a  pK  of  4.0,  with  the  fluoride  ion  binding  faster  to 
the  protonated  species. 

It  was  felt  that  a  similar  study  of  the  binding  of 
cyanide  to  HRP  over  a  wide  pH  range  might  yield  additional 
information  about  the  nature  of  ligand  binding  to  HRP. 

Hence,  after  the  construction  of  the  apparatus  described  in 
the  next  chapter,  such  a  project  was  initiated. 


, 
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CHAPTER  2. 

APPARATUS  FOR  THE  STUDY  OF  FAST  REACTIONS  IN  SOLUTION 

•»  f, 

The  title  of  this  chapter  raises  a  questions  how  fast 
is  a  ”fast  reaction”?  Caldin43  has  defined  a  fast  reaction 
as  ” • . «  a  reaction  that  is  fast  relative  to  the  time  re¬ 
quired  for  mixing  and  observation  by  conventional  methods.” 
By  Caldin’s  own  admission,  this  definition  is  imprecise,  but 
serviceable.  If  one  considers  the  conventional  method  of 
mixing  reagents  in  a  cuvet  and  taking  a  reading  on  a  spec¬ 
trophotometer,  a  fast  reaction  would  then  be  one  with  a 
half-time  of  less  than  a  few  seconds.  To  study  such  a 
reaction,  one  needs  a  special  apparatus.  By  constructing 
an  apparatus  designed  to  mix  reagents  rapidly,  one  can  study 
reactions  with  a  minimum  half-time  of  a  few  milliseconds. 

For  even  faster  reactions,  one  must  use  a  method  which 
avoids  mixing  altogether.  Such  a  method  is  the  temperature- 
jump  technique,  which  is  one  of  the  relaxation  methods  pro¬ 
posed  by  Eigen43,  all  of  which  involve  taking  a  reversible 
reaction  system  which  is  at  equilibrium  and  suddenly  per¬ 
turbing  the  equilibrium  so  that  the  reaction  will  then 
proceed  until  a  new  position  of  equilibrium  is  established. 
The  time  taken  to  reach  the  new  position  of  equilibrium  can 
be  related  to  the  rate  constants  of  the  forward  and  reverse 
reactions.  The  use  of  relaxation  techniques  makes  it  possi¬ 
ble  to  study  reactions  with  half-times  as  small  as  a  few 
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microseconds  . 

This  chapter  contains  a  discussion  of  the  construction 
and  performance  of  two  types  of  apparatus  for  the  study  of 
fast  reactions  in  solution:  a  stopped-flow  apparatus  and  a 
temperature- jump  apparatus. 

Stopped-Flow  Apparatus 
History  and  Principles 

Rapid  mixing  and  flow  techniques  date  back  to  1923 
when  Hartridge  and  Roughton  constructed  an  apparatus  which 
achieved  mixing  times  of  one  millisecond44.  Because  of  the 
slow  response  time  of  the  detection  instrumentation  avail¬ 
able  then,  time  resolution  was  achieved  by  a  continuous  flow 
procedure  in  which  the  freshly  mixed  solution  passed  down  an 
observation  tube  and  the  composition  of  the  reaction  mixture 
was  determined  at  various  cross  sections  by  optical,  ther¬ 
mal,  or  electrical  analysis.  If  the  flow  velocity  is  kept 
constant,  the  distance  from  the  mixing  chamber  will  corre¬ 
spond  to  a  time  equal  to  the  distance  divided  by  the  ve¬ 
locity. 

The  stopped-flow  method  was  developed  by  Chance  in 
19404 6,  but  has  only  come  into  widespread  use  since  Gibson’s 
modifications  in  195246.  This  method  involves  rapid  mixing 
and  flow  of  mixed  reactants  to  an  observation  chamber.  The 
flow  is  stopped  quickly  and  the  progress  of  the  reaction  in 
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the  observation  chamber  is  followed  by  some  appropriate 
method.  The  stopped- flow  technique  has  several  advantages 
over  the  continuous  flow  method,  including  the  fact  that 
the  observation  is  independent  of  the  rate  and  character  of 
flow  down  the  observation  tube.  Also,  very  small  volumes 
of  reactants  can  be  used,  and  a  permanent  record  of  the 

*  • 

progress  of  the  reaction  throughout  its  entire  course  can 
be  obtained. 

Description  of  the  Apparatus 

Ideas  for  the  design  of  the  stopped-flow  apparatus  were 
obtained  from  the  papers  of  Gibson47  and  of  Sturtevant48 
and  from  examination  of  Dr.  R.  A.  Alberty’s  stopped-flow 
at  the  University  of  Wisconsin,  which  was  built  by  D.  A. 

Ver  Ploeg49.  The  primary  innovation  introduced  in  the 
model  described  in  this  chapter  involves  a  new  design  for 
the  mixing  chamber  and  observation  section.  A  block  dia¬ 
gram  of  the  apparatus  is  shown  in  Figure  4. 

A  spectrophotometric  detection  system  is  used  which 
consists  of  a  Bolex  Ace  50  watt  projection  lamp;  a  Bausch 
and  Lomb  No.  33-86-26-07  High  Intensity  Model  5  grating 
monochromator  with  a  Bausch  and  Lomb  quartz  collective  lens 
No.  1;  a  cylindrical  lens  made  from  a  Fisher  plastic  ther¬ 
mometer  magnifier;  an  R.C.A.  Model  1P28  photomultiplier  in 
a  circuit  with  an  emitter  follower  output;  and  a  Tektronix 
535^  oscilloscope  with  a  Type  D  high  gain  differential 
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Figure  4.  Block  diagram  of  stopped- flow  apparatus 
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plug- in  preamplifier.  A  Hewlett-Packard  Harrison  Model 

6264a  dc  power  supply  for  the  tungsten  lamp  provides  eight 
volts  with  0.01  °/0  line  and  load  regulation  and  500  juV 
ripple.  Photomultiplier  voltage  is  obtained  from  a  Hewlett- 
Packard  Harrison  Model  6110A  dc  power  supply.  Experimental 
traces  are  photographed  with  a  35  mm  camera  mounted  on  the 
oscilloscope  with  a  Beattie-Coleman  Polexa  adapter. 

The  photomultiplier  circuit  is  shown  in  detail  in 
Figure  5  and  its  power  supply  unit  in  Figure  6.  These  two 
circuits,  built  by  the  University  of  Alberta  Chemistry 
Electronics  Shop,  contain  a  number  of  useful  features.  The 
sensitivity  can  be  varied  by  selecting  a  value  of  load  re¬ 
sistance  ranging  from  30  k-^  to  10  M/x.  In  order  to  filter 
the  noise  from  the  output  signal,  an  RC  circuit  of  the  ap¬ 
propriate  time  constant  can  be  introduced  by  combining  the 
load  resistor  with  a  variable  capacitance  ranging  from  0  to 
1  pF.  The  circuit  also  contains  a  zero  offset  control  to 
provide  various  amounts  of  bucking  voltage  to  the  photo¬ 
multiplier  output.  This  allows  the  initial  photomultiplier 
signal  to  be  nulled  and  thus  provides  a  convenient  means  for 
observing  a  small  signal  change  which  is  imposed  upon  a 
larger  signal. 

In  the  case  of  an  experiment  involving  a  reaction 
which  proceeds  slowly  enough  that  an  oscilloscope  sweep  rate 
of  slower  than  0.1  sec/cm  can  be  used,  the  procedure  for  ob¬ 
taining  a  record  of  the  oscilloscope  trace  is  very  simple. 
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Figure  5.  Photomultiplier  circuit 
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*These  resistors  are  selected  so  that  the  meter  reads  in 
"Batt.  O.K.”  sector  (0.6-0. 8  ma )  while  batteries  are 
under  load. 


Figure  6.  Power  supply  for  photomultiplier  circuit 
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The  appropriate  sweep  rate  is  chosen  and  the  sweep  is  set 
for  free-running  operation.  When  the  trace  appears  on  the 
left  side  of  the  oscilloscope  screen,  the  camera  shutter  is 
opened  and  the  driving  piston  of  the  stopped- flow  apparatus 
is  pushed  in  quickly.  When  the  reaction  reaches  equilibrium, 
as  evidenced  by  no  further  change  in  trace  voltage  with 
time,  the  camera  shutter  is  closed.  For  faster  reactions, 
a  trigger  is  necessary.  As  shown  in  Figure  4,  the  trigger 
consists  of  a  45  V  battery  in  a  circuit  which  is  completed 
when  the  driving  piston  hits  the  stop.  Generally  it  is 
important  to  observe  the  first  portion  of  the  trace  and  to 
be  able  to  establish  unequivocally  the  point  corresponding 
to  zero  reaction  time.  These  aims  are  best  achieved  if  a 
small  wire  is  attached  to  the  stop  in  such  a  manner  that 
the  piston  face  brushes  the  wire  just  before  hitting  the 
stop.  Thus,  the  oscilloscope  is  triggered  at  such  a  time 
that  the  last  portion  of  the  flow,  and  hence,  the  first  part 
of  the  reaction  are  always  observed. 

The  overall  spectral  sensitivity  of  the  detection  sys¬ 
tem  is  a  function  of  the  relative  spectral  composition  of 
the  lamp,  the  transmittance  curve  of  the  monochromator 
grating,  and  the  spectral  response  of  the  photomultiplier. 

It  is  important  to  know  the  spectral  sensitivity  of  the 
system,  since  a  reaction  is  often  encountered  which  can  be 
followed  at  more  than  one  wavelength  of  light.  If  one  se¬ 
lects  a  wavelength  corresponding  to  a  high  sensitivity, 
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then  a  lower  photomultiplier  voltage  can  be  used  to  get  the 
same  output  voltage,  and  a  reduction  of  photomultiplier 
voltage  results  in  a  decrease  of  shot  noise  and  hence  in  a 
better  signal- to-noise  ratio*  Such  a  sensitivity  curve  for 
the  detection  system  described  above  is  shown  in  Figure  7. 

It  was  obtained  by  filling  the  stopped- flow  observation  tube 
with  water  and  measuring  the  output  voltage  of  the  photo¬ 
multiplier  vs.  wavelength*  A  calibrated  oscilloscope  was 
used  for  the  voltage  measurements* 

The  stopped-flow  device  itself  is  best  described  in 
conjunction  with  a  photograph  of  the  apparatus,  Figure  8. 

The  main  body  measures  10*5  cm  x  43  cm  and  is  made  of 
Plexiglas,  painted  black  to  keep  out  stray  light.  The 
driving  piston  (a)  is  made  of  brass,  and  it  slides  through 
two  Plexiglas  collars*  The  stopping  pin  (b)  is  made  of 
stainless  steel  and  it  seats  1*5  cm  into  the  main  body  of 
the  apparatus.  There  are  five  positions  for  the  stop, 
spaced  1*05  cm  apart*  The  thermostated  compartment  (C ) , 
shown  with  the  top  off,  is  4*0  cm  deep*  It  contains  the 
two  1*0  ml  driving  syringes  (D)  which  lead  to  Hamilton 
three-way  valves  (E)*  The  valves  are  positioned  so  that 
reactants  can  flow  into  the  driving  syringes  from  thermo¬ 
stated  5*0  ml  reservoir  syringes  (F)  or  from  the  driving 
syringes  through  Teflon  tubing  to  the  mixing  and  observa¬ 
tion  chamber  (G).  From  there  the  used  reaction  mixture 
passes  to  the  2.0  ml  collection,  or  back-up  syringe  (H). 


- 
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Figure  7 • 


Spectral  sensitivity  of  optical  detection  system 
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Figure  8.  Photograph  of  Stopped- Flow  Apparatus 

Key:  A,  driving  piston;  B,  stopping  pin;  C, 
thermostated  compartment;  D,  driving  syringes; 
E,  three-way  valves;  F,  reservoir  syringes; 

G,  mixing  and  observation  chamber;  H,  back-up 
syringe;  I,  monochromator;  J,  photomultiplier. 
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This  syringe  provides  a  small  back  pressure  to  help  overcome 
the  inertia  of  the  driving  syringes  .  Another  Hamilton 
three-way  valve  allows  the  back-up  syringe  to  be  emptied 
easily,  or  bypassed  if  desired.  Also  shown  in  the  photo¬ 
graph  are  the  monochromator  (i)  and  the  box  housing  the 
photomultiplier  circuit  (j). 

The  heart  of  the  stopped- flow  apparatus  is  the  mixing 
and  observation  section.  Figure  9^,  drawn  to  scale,  shows 
a  top  exploded  view  of  this  section,  which  consists  of  a 
series  of  wafers  and  light  ports,  all  made  of  Plexiglas. 

The  arrows  show  the  paths  of  the  reactants.  The  two  halves 
of  this  section  are  identical  except  that  the  effluent  tube 
is  located  in  one  of  the  nBn  wafers  rather  than  the  center. 
The  light  ports  are  hollow  Plexiglas  rods  with  a  0.95  cm 
o.d.  and  a  0.55  cm  i.d.  They  are  glued  to  wafer  A  on  one 
side  and  on  the  other  side  attach  to  the  thermostated  com¬ 
partment  by  means  of  Plexiglas  nuts.  Wafer  A  contains  a 
viewing  window  made  from  a  rod  of  General  Electric  type  104 
fused  quartz,  5  nun  in  diameter  and  2.0  cm  in  length.  The 
ends  of  the  rod  are  ground  and  polished  to  a  plate  glass 
finish.  The  window  is  held  in  place  by  epoxy  resin.  Be¬ 
cause  the  Teflon  tubing  cannot  be  glued  in  place,  a  special 
connection  is  necessary;  it  is  designed  so  that  the  flared 
tubing  fits  over  a  stainless  steel  fitting  which  is  glued 
into  the  Plexiglas.  The  tubing  is  held  snugly  in  place  by 
a  hollow  threaded  steel  rod.  The  length  of  the  observation 
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Figure  9»  Stopped-flow  mixing  and  observation  section. 
A.  Exploded  top  view,  to  scale.  B.  Outside  face  of 
wafer  B.  C.  Inside  face  of  wafer  B.  B  and  C  are  en¬ 


larged  by  a  factor  of  two 
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tube,  jl.  e.  two  times  the  width  of  a  B  wafer,  is  loO 67  cm. 

The  two  faces  of  wafer  B  are  illustrated,  enlarged  by 
a  factor  of  two,  in  Figures  9B  and  9C.  Solution  containing 
one  of  the  reactants  flows  into  the  center  of  the  splitter 
shown  in  the  drawing  of  the  outside  face.  The  solution  is 
then  split  into  four  portions  which  pass  through  the  wafer 
and  meet  head-on  with  the  streams  of  the  other  reactant  so¬ 
lution.  This  constitutes  the  first  stage  of  mixing.  The 
four  freshly  mixed  portions  flow  down  the  inside  face  of  the 
wafer  and  mix  a  second  time  when  they  meet  tangentially  at 
the  center  cavity.  From  there  the  solution  passes  through 
a  channel  to  the  center  of  the  observation  tube,  where  it 
splits  to  fill  the  observation  tube.  The  line  of  flow  con¬ 
tinues  out  of  the  observation  tube  to  the  back  of  the  mix¬ 
ing-observation  section  and  then  into  a  stainless  steel 
hypodermic  needle  which  connects  to  the  three-way  valve  at 
the  back  of  the  apparatus. 

The  wafers  are  held  together  by  four  screws,  one  placed 
above  and  one  below  the  tubing  inlets  and  the  other  two  at 
corresponding  positions  at  the  rear  of  the  section.  Under 
high  pressure  driving,  a  few  small  leaks  appeared  from  be¬ 
tween  the  wafers.  These  were  stopped  by  coating  the  whole 
mixing-observation  section  with  General  Electric  silicone 
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Performance  of  the  Apparatus 

The  first  performance  criterion  a  stopped-f low  appara¬ 
tus  must  meet  is  that  it  must  have  the  ability  to  flush  old 
solution  completely  from  the  observation  chamber.  For  the 
apparatus  described  above,  the  total  volume  of  the  mixing 
chamber  and  observation  tube  is  about  0.07  nil  and  the  volume 
delivered  by  the  movement  of  the  driving  piston  from  one 
stop  to  the  next  is  0.4  ml.  Thus,  since  at  least  four  vol¬ 
umes  of  fresh  solution  sweep  out  the  observation  tube  before 
the  portion  which  will  be  observed  enters,  it  is  expected 
that  only  a  negligible  amount  of  old  solution  would  remain 
behind.  This  expectation  is  borne  out  by  experiments  using 
colored  solutions  and  photometric  detection  which  showed 
that  no  detectable  old  solution  remained  in  the  observation 
chamber  when  a  portion  of  fresh  solution  was  injected. 

Another  performance  criterion  is  the  time  resolution, 

i 

which  depends  on  two  factors  *  the  mixing  time  and  the  time 
for  flow  from  the  mixing  chamber  to  the  observation  tube, 
the  latter  being  the  ’’dead  time”.  The  mixing  time  is  usual¬ 
ly  fast  compared  to  the  dead  time,  and  this  can  be  confirmed 
by  performing  a  simple  experiment.  Since  an  acid-base  reac¬ 
tion  is  instantaneous  on  the  stopped- flow  time  scale,  any 
change  in  transmittance  with  time  after  an  acid  and  base  are 
mixed  in  the  presence  of  an  indicator  must  be  due  to  incom¬ 
plete  mixing  at  the  beginning.  For  this  test,  one  syringe 
was  filled  with  a  dilute  solution  of  brom  cresol  green 


' 
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indicator  >±n  base  form,  and  the  other  with  a  solution  of  pH 
2.0  buffer .  The  indicator  changes  from  blue  to  yellow  as 
the  solution  becomes  acidic  through  the  pH  range  3 • 8-5*4. 

The  monochromator  was  set  at  620  mp,  which  is  the  wavelength 
corresponding  to  the  maximum  of  the  absorption  peak  of  the 
base  form  of  the  indicator,  the  reactants  mixed  and  the  re¬ 
sulting  oscilloscope  trace  was  photographed.  Inspection  of 
the  photograph  showed  that  the  solution  was  greater  than 
98  °/ o  mixed  by  the  time  the  solution  was  stopped  in  the 
observation  chamber.  Thus,  the  time  resolution  was  shown  to 
depend  primarily  on  the  dead  time  of  the  apparatus. 

The  dead  time  of  the  apparatus  can  be  calculated  if  the 
flow  rate  and  volume  of  the  mixing  chamber  and  observation 
tube  are  known  exactly.  If  these  are  not  known,  the  dead 
time  can  be  found  easily  by  experiment  .  The  point  on  the 
experimental  oscilloscope  trace  corresponding  to  the  time 
that  the  flow  is  stopped  is  the  apparent  t=0.  The  true  t=0 
exists  when  the  reactants  are  mixed,  so  it  is  the  time  be¬ 
tween  the  true  and  apparent  t=0  which  must  be  determined. 

The  reaction  studied  to  determine  this  time  was  the  redox 
reaction  between  ascorbate  and  f erricyanide.  It  was  fol¬ 
lowed  at  422  mp,  the  maximum  of  the  absorption  peak  of 

The  author  is  indebted  to  members  of  Dr.  G.  G.  Hammes 1 
group  at  Cornell  University  for  helpful  discussions  regard¬ 
ing  this  experiment. 
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ferricyanide.  A  solution  of  2xlO"4M  potassium  ferricyanido 
was  mixed  with  excess  ascorbate  and  photographs  were  taken 
of  the  resulting  oscilloscope  traces.  Three  experiments 
were  performed  using  the  same  ferricyanide  concentration, 
but  each  with  a  different  ascorbate  concentration.  As  long 
as  the  ascorbate  is  kept  in  excess  so  that  the  reaction  is 
pseudo-first- order ,  it  is  not  necessary  to  know  the  exact 
ascorbate  concentrations  to  carry  out  the  analysis. 

Photographs  of  the  oscilloscope  traces  provided  a  record 
of  photomultiplier  circuit  output  voltage  vs.  time.  This 
output  voltage  is  proportional  to  the  current  at  the  anode 
of  the  photomultiplier  tube,  which  is  in  turn  proportional 
to  the  intensity  of  light  incident  upon  the  photocathode. 

If  I0  and  I  are  the  light  intensities  with  solvent  only  and 
with  solvent  plus  absorbing  compound,  respectively,  and  V0 
and  V  are  the  corresponding  voltages  which  are  produced, 
thens  V/Vo  =  i/lo;  but  i/lo  =  T,  the  transmittance,  so  V/V0 
is  also  equal  to  the  transmittance.  To  analyze  the  results 
of  the  dead  time  determination  experiments,  the  oscillo¬ 
scope  screen  was  calibrated  in  terms  of  transmittance  using 
solutions  whose  transmittance s  were  known  from  measurement  on 
a  Beckman  DU  spectrophotometer.  Points  were  read  from  the 
photographs  with  the  aid  of  a  photographic  enlarger,  con¬ 
verted  to  transmittance  and  then  absorbance  readings » 

For  a  first-order  reaction: 

kt 

2.303 


-log  c  =  -log  Co  + 


(1) 
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where  Co  and  c  are  the  initial  concentration  and  the  concen¬ 
tration  at  time  t  of  the  first-order  reactant.  Since  by 
Beer’s  law,  A<*c,  a  plot  of  -log  A  vs.  t  should  yield  a 
straight  line  with  an  intercept  at  -log  Ao»  Figure  10  con¬ 
tains  such  plots  for  the  dead  time  test,  with  -log  A  of  fer- 
ricyanide  plotted  vs.  apparent  time  for  the  three  different 
concentrations  of  ascorbate.  If  no  reaction  had  taken  place 
at  the  apparent  t=0,  the  three  lines  should  have  intersected 
at  that  point,  since  the  true  c0  of  ferricyanide  and  hence, 
the  true  log  Ao,  would  be  the  same  for  the  three  experiments. 
This  was  not  observed,  as  shown  in  Figure  10.  Extrapolation 
of  the  three  lines  back  to  the  point  where  they  intersect, 
and  thus  at  which  they  do  represent  equal  values  of  -log  A0, 
allows  determination  of  the  true  t=0.  From  Figure  10,  it 
can  be  seen  that  the  true  t=0  occurs  approximately  ten  msec 
before  the  apparent  t=0,  so  the  dead  time  of  the  apparatus 
is  about  ten  msec.  The  small  locus  of  intersection  indi¬ 
cates  that  although  the  driving  piston  was  pushed  by  hand, 
the  flow  rate  achieved  was  quite  reproducible.  It  is  ex- 

i 

pected  that  a  shorter  dead  time  would  be  possible  if  the 
driving  piston  were  driven  pneumatically  to  increase  the 
flow  rate. 


Temperature- Jump  Apparatus 
History  and  Principles 


The  general  theory  of  relaxation  methods  is  treated  in 
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Apparent  Time  (msec) 


Figure  10.  Plot  of  -log  A  of  ferricyanide  vs.  appar¬ 
ent  reaction  time.  The  negative  value  of  time  at 

t 

which  the  lines  intersect  is  the  dead  time  of  the 
stopped-flow  apparatus. 
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detail  elsewhere60"53,  so  this  section  will  deal  only  with 
the  temperature- jump  method  and  apparatus. 

The  temperature- jump  method  was  first  proposed  by  Eigen 


in  195443 ,  and  a  description  of  the  first  operational  instru¬ 


ment  appeared  five  years  later53.  Since  that  time,  the  use 
of  this  technique  has  grown  very  rapidly,  and  has  found  wide¬ 
spread  use  in  the  study  of  inorganic  and  biochemical  reac¬ 
tions. 

The  thermodynamic  basis  for  the  temperature- jump  method 
is  expressed  by  a  form  of  the  Gibbs-Helmholtz  equation: 


AH° 

RT3 


(2) 


From  eq  2  it  can  be  seen  that  the  constant  pressure  equilib¬ 
rium  constant  of  a  reaction  will : change  upon  a  perturbation 
of  the  temperature,  provided  the  reaction  being  studied  has 
a  non-zero  standard  enthalpy.  A  reaction  which  has  a  very 
small  value  of  AH0  can  often  be  studied  by  coupling  it  to 
one  which  has  a  greater  value  of  ah0. 

The  temperature  rise  is  obtained  by  discharging  a  capac¬ 
itor  through  a  solution  containing  the  reactants  at  equilib¬ 
rium.  There  is  no  spark  or  flash  through  the  solution;  the 
rise  in  temperature  is  due  to  a  Joule  heating  of  the  solu¬ 
tion.  This  heating  may  be  likened  to  the  heating  caused  by 
an  IR  drop  across  a  resistor.  In  the  solution,  the  energy 
is  distributed  among  the  translational,  vibrational,  and 
rotational  degrees  of  freedom  in  a  time  which  is  much  faster 
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than  the  time  required  for  the  discharge  of  the  capacitor64. 
Since  ions  must  be  present  to  carry  the  current,  the  solu¬ 
tion  is  usually  made  0.1  M  in  some  inert  electrolyte  such  as 
potassium  nitrate. 

The  magnitude  of  the  temperature  rise,  AT,  can  be  esti¬ 
mated  by  assuming  that  all  the  energy  stored  in  the  capaci¬ 
tor  is  transferred  to  the  solution,  and  is  given  by: 


AT  = 


CV* 

2p  vk 


(3) 


where  C  is  the  capacitance  in  farads,  V  is  the  voltage  on 
the  capacitor  in  volts >  p  is  the  solution  density  in  g/cm3, 
v  is  the  volume  of  liquid  heated,  in  cm3,  and  k  is  the  spe¬ 
cific  heat  capacity  of  the  solution  in  J/g  deg.  A  typical 
temperature- jump  apparatus  is  designed  to  produce  tempera¬ 
ture  increases  of  two  to  ten  degrees. 

The  temperature- jump  method  is  applicable  over  a  wide 
time  scale.  The  lower  limit  is  determined  by  the  time  re¬ 
quired  to  heat  the  solution.  This  is  usually  expressed  in 
terms  of  a  heating  time  constant,  Tc»  If  the  inductance  of 
the  capacitor  and  the  circuit  capacitance,  inductance,  and 
resistance  are  assumed  to  be  negligible,  T 0  may  be  calcu¬ 
lated  from  the  equation: 


To  =  Va  RC 


(4) 


where  R  is  the  resistance  of  the  solution  and  C  is  the  capac¬ 
itance  of  the  capacitor.  The  values  of  To  for  various  tem¬ 
perature-jump  instruments  range  from  0.1  p. sec  to  10  /asec66. 
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The  upper  time  limit  for  the  apparatus  is  of  the  order  of 
one  second.  This  limit  is  associated  with  the  cooling  of 
the  solution  between  the  electrodes  by  convection  currents 
from  the  bulk  solution. 

Description  of  the  Apparatus 

The  temperature- jump  apparatus  described  here  is  mod¬ 
eled  after  the  one  used  by  Dr.  R.  A.  Alberty  and  co-workers 
at  the  University  of  Wisconsin66. 

The  detection  system  used  with  the  instrument  is  like 
the  one  illustrated  in  Figure  4  and  described  in  the  pre- 
ceeding  section,  but  with  the  cylindrical  lens  replaced  by  a 
spherical  lens  from  a  Parr  Instrument  Co,  buret  meniscus 
magnifier.  This  lens  was  chosen  because  it  fit  conveniently 
into  the  light  port  of  the  temperature- jump  cell  thermo- 
stating  jacket. 

A  schematic  diagram  of  the  temperature- jump  circuit  is 
shown  in  Figure  11.  A  Plastic  Capacitors,  Inc.  50  kV,  5 
High  Voltage  Power  Supply  Model  HV-500-502M  coupled  with  a 
Superior  Power stat  No.  10B  is  used  to  charge  a  Plastic  Capac¬ 
itors,  Inc.  LK500-104ZMD,  0.1  pF,  50  kV  capacitor.  The  ca¬ 
pacitor  is  discharged  through  an  EG  and  G  Inc.  Spark  Gap 
Model  GP14B  by  triggering  the  spark  gap  with  an  EG  and  G  Inc. 
Trigger  Module  Model  TM-11.  The  oscilloscope  is  triggered 
by  means  of  a  signal  inductively  produced  in  a  coil  around 
the  lead  from  the  capacitor.  A  second  temperature- jump 
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Figure  11 


Schematic  diagram  of  temperature- jump  system 
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apparatus  was  built  which  differed  from  the  one  described 
above  in  that  the  spark  gap  was  constructed  by  the  Univer¬ 
sity  of  Alberta  Chemistry  Electronics  Shop  and  the  high 
voltage  is  supplied  by  a  Spellman  High  Voltage  Power  Supply 
Model  2040. 

The  features  of  the  temperature- jump  cell  are  illus¬ 
trated  in  Figure  12,  a  side  sectional  view  of  the  cylin¬ 
drical  cell,  drawn  to  scale.  The  main  body  of  the  cell  is 
made  of  Plexiglas,  the  base  of  Colorlith,  the  windows  of 
quartz  rod  as  described  in  the  preceeding  section,  and  the 
top  and  electrodes  of  stainless  steel.  All  corners  of  the 
electrodes  were  rounded  and  the  surfaces  polished  to  prevent 
arcing.  The  volume  of  solution  contained  between  the  elec¬ 
trode  faces  is  1  cm3,  and  experiments  can  be  performed  us¬ 
ing  as  little  as  5  ml  total  solution  in  the  cell. 

The  cell  is  fitted  snugly  into  a  brass  thermostating 
jacket?  the  cell  and  jacket  are  both  surrounded  by  a  light¬ 
tight  canister.  The  high  voltage  connection  is  made  by  the 
banana  plug  at  the  bottom  of  the  cell,  and  the  top  electrode 
is  connected  to  ground. 

Performance  of  the  Apparatus 

The  resistance  of  the  cell  containing  0.11  M  KNO3 , 
measured  with  a  bridge  circuit,  is  133  ohms.  Using  this 
value  in  eq  4,  along  with  a  capacitance  of  0.1  pF,  the  heat¬ 
ing  time  constant  for  the  apparatus  is  calculated  to  be 
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Temperature- jump  cell 
stainless  steel 
sty  plexiglas 
quartz 
colorlith 


Figure  12. 
Key : 
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7psec.  Although  this  figure  represents  the  minimum  relaxa¬ 
tion  time  which  can  be  studied  using  this  apparatus,  one 
must  also  consider  the  response  time  of  the  photomultiplier 
circuit.  The  response  time  of  the  circuit  depends  on  the 
values  of  the  residual  circuit  capacitance,  the  filtering 
capacitance,  and  the  load  resistance.  Values  of  the  circuit 
response  time  constant,  7^c»  were  determined  experimentally. 

A  xenon  flash  lamp  with  a  flash  decay  time  constant  of  3 
psec  was  placed  in  front  of  the  photomultiplier  tube.  With 
the  filtering  capacitance  set  at  zero,  photographs  of  oscil¬ 
loscope  traces  of  the  photomultiplier  circuit  response  fol¬ 
lowing  a  lamp  flash  were  taken  at  each  setting  of  the  load 
resistor.  7^c  was  obtained  from  inspection  of  the  photo¬ 
graphs,  and  was  taken  to  be  the  time  required  for  the  output 
voltage  to  decay  to  two-thirds  of  the  maximum  value  caused 
by  the  lamp  pulse.  Table  II  contains  values  of  7^c  for  the 
various  values  of  load  resistance  in  the  photomultiplier 
circuit. 


Table  II 

Response  Time  Constants  for  the  Photomultiplier  Circuit 

.  _  .  .  n- 


Load  Resistance  / c.  .  jusec 

30  k^  3 

100  k  jy-  5 

300  k-"-  12 

1  Ma  30 

3  M-^  60 

10  125 
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The  theoretical  expression  for  the  temperature  rise  of 
the  solution,  given  by  eq  3,  can  be  simplified  for  the  usual 
conditions  of  operation  for  this  apparatus.  With  v  equal  to 
1  cm3,  C  equal  to  0.1  pF,  and  values  of  k  and  j>  taken  as 
4.18  and  1.0,  the  corresponding  values  for  water,  eq  3  be¬ 
comes  : 

AT  =  1.19  x  lO*”8  V3  (5) 

The  validity  of  eq  5  was  checked  by  another  worker  in  this 
laboratory67  using  buffers  of  known  AH°  and  a  pH  indicator 
according  to  the  method  described  by  Czerlinski 68 .  Experi¬ 
mental  values  of  A  T  agreed  within  5  °/o  of  the  values  pre¬ 
dicted  by  eq  5 • 


Combination  Apparatus  Table 

Because  both  the  stopped-flow  and  temperature- jump 
instruments  use  the  same  detection  system,  it  was  felt  that 
it  would  be  convenient  to  be  able  to  use  the  two  instruments 
interchangeably.  Two  identical  tables  were  built  for  this 
purpose.  Figure  13  is  a  photograph  of  one  of  the  tables, 
set  up  for  the  stopped-flow  apparatus,  labeled  nC"  on  the 
photograph.  The  frame  of  the  table  is  constructed  of  2" 
angle  iron,  with  additional  support  from  the  electronic  com¬ 
ponent  cabinet  on  the  left  side.  The  top  of  the  table  and 
the  shelf  on  the  right  bottom  are  made  of  3/k"  Colorlith. 

The  lamp  (A) ,  monochromator  (b)  and  photomultiplier  (D)  are 
held  by  clamps  which  slide  along  an  optical  bench  (E)  made 
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Figure  13 •  Apparatus  Table 

Key:  A,  lamp;  B,  monochromator;  C,  stopped- 
flow  apparatus;  D,  photomultiplier;  E,  optical 
bench;  F,  oscilloscope;  G,  camera  and  adapter; 
H,  constant  temperature  bath;  I,  box  containing 
capacitor  and  spark  gap;  J,  trigger  output; 

K,  high  voltage  discharge  switch;  L,  trigger 
module.  The  cabinet  on  the  left  contains  the 


various  power  supplies. 
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of  two  pieces  of  ln  angle  iron  spaced  3/^n  apart .  Details 
of  the  optical  bench  and  mounting  arrangement  can  be  seen  in 
more  detail  in  Figure  8.  Also  shown  in  Figure  13  are  the 
oscilloscope  (f),  camera  and  adapter  (g)  and  constant  tem¬ 
perature  bath  (h).  The  capacitor  and  spark  gap  are  con¬ 
tained  in  the  metal  box  (i),  located  behind  the  trigger  out¬ 
put  (J),  high  voltage  discharge  switch  (K)  and  trigger  mod¬ 
ule  (L).  The  cabinet  contains  the  various  power  supplies 
and  the  batteries  for  the  photomultiplier  zero  offset  con¬ 
trol.  Conversion  from  a  stopped-flow  apparatus  to  a  tem¬ 
perature-jump  apparatus  is  simply  a  matter  of  unbolting  the 
one,  replacing  it  with  the  canister  housing  the  temperature- 
jump  cell,  and  checking  the  alignment  of  the  optical  system, 
a  procedure  which  takes  only  a  few  minutes. 
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CHAPTER  3.  THE  KINETICS  OF  CYANIDE  BINDING 

BY  FERRIC  HORSERADISH  PEROXIDASE 

A  study  of  the  kinetics  of  cyanide  binding  to  HRP 
might  be  expected  to  yield  fundamental  information  about 
the  nature  of  the  active  site  in  peroxidase.  An  early  study 
of  this  reaction  at  pH  4.2  and  6.2  showed  little  change  in 
the  value  of  the  forward  rate  constant38.  This  chapter  con¬ 
tains  a  report  of  a  study  of  this  reaction  at  nineteen  val¬ 
ues  of  pH 9  ranging  from  4.21  to  11.31* 

Experimental 

Kinetic  measurements  were  made  using  the  stopped- flow 
apparatus  and  spectrophotometric  detection  system  described 

£- 

Abbreviations  used  in  this  chapter:  HRP  or  P,  ferric 
horseradish  peroxidase;  RZ,  purity  number,  after  Theorell 
and  Maehly14;  p,  ionic  strength;  CN  and  PCN,  all  ionizable 
species  of  cyanide  and  complex;  HP  and  HgP,  protonated  forms 
of  peroxidase;  [P],  [CliNfJ,  and  [PCN],  equilibrium  concentra¬ 
tions  of  all  forms  of  peroxidase,  cyanide,  and  complex;  [P]o 

and  [CN]o j  initial  concentrations  of  all  forms  of  peroxidase 
and  cyanide;  KL,  dissociation  constant  of  ligand;  K^,  disso¬ 
ciation  constant  of  peroxidase-cyanide  complex;  ki  and  kli, 
apparent  second-order  binding  rate  constant  and  apparent 
first-order  dissociation  rate  constant. 
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in  Chapter  2.  Auxiliary  absorbance  measurements  and  spectro- 
photometric  equilibrium  constant  determinations  were  made  on 
a  Beckman  DU  spectrophotometer,,  The  cell  compartment  of  the 
DU  and  the  stopped- flow  apparatus  were  both  thermostated  at 
25 » 0*0.1  °C.  A  Beckman  expanded  scale  pH  meter  was  used  for 
all  pH  measurements. 

HRP  (analytical  reagent  grade,  RZ  =  2„9j  lot  648520 6 
and  purified  grade,  RZ  =  0.9 )  was  obtained  from  the 
Boehringer -Mannheim  Corp. ,  N.  Y.  Reagent  grade  potassium 
cyanide  (Fisher)  was  used  without  further  purification.  The 
concentrations  of  stock  KCN  solutions  were  determined  by 
titration  with  a  standard  silver  nitrate  solution  using 
potentiometric  determination,  of  the  end  point.  Below  pH  8, 
possible  loss  of  HCN  was  minimized  by  using  solutions  of 
cyanide  which  were  less  than  10“ 4  M  and  by  keeping  them  in 
tightly  closed  vessels  which  contained  only  a  small  volume 
of  airc  All  solutions  were  made  from  conductivity  water. 

Solutions  of  HRP  and  cyanide  were  each  made  to  the  de¬ 
sired  pH  using  buffer  of  ionic  strength  equal  to  0.01  and 
sufficient  potassium  nitrate  so  that  the  total  ionic  strength 
was  constant  at  p.  equal  to  0.11.  Concentrations  of  HRP  and 
cyanide  were  varied  at  different  pH  values  in  such  a  manner 
that  the  experimental  absorbances  ranged  from  0.25  to  0.50, 
which  is  in  the  region  of  minimum  error  due  to  instrumental 
uncertainty59.  The  total  HRP  concentration  was  typically 
5  x  10~6  M,  while  the  total  cyanide  concentration  ranged 
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from  about  5  x  10" 6  M  to  1  x  10" 4  M  over  the  pH  range  in 
which  experiments  were  performed. 

The  buffer  systems  employed  at  various  pH  ranges  were: 
acetate  below  pH  5»5»  cacodyl ic  acid  between  pH  5*5  and  7*0, 
tris-nitric  acid  between  pH  7.0  and  9«0»  borate  between  pH 
8.0  and  10.0,  glycine  between  pH  9°0  and  11.0,  and  potassium 
hydroxide  above  pH  11.  Each  buffer  was  checked  for  possible 
complexing  with  HRP  by  examining  the  spectrum  of  HRP  at  a 
particular  pH  in  the  presence  and  absence  of  buffer.  In 
addition,  experiments  were  repeated  with  two  or  more  differ¬ 
ent  buffers  in  all  overlapping  regions.  In  all  cases,  ex¬ 
periments  at  the  same  pH  but  with  different  buffers  gave 
results  agreeing  within  experimental  error. 

HRP  and  its  complexes  are  known  to  undergo  splitting 
at  extremes  of  pH30.  Experiments  at  the  pH  extremes  were 
performed  with  sufficient  rapidity  that  the  decrease  in 
absorbance  of  HRP  at  403  mju  during  the  experiment  was  less 
than  1  °/0o  Above  pH  10  problems  in  pH  control  were  encoun¬ 
tered  because  of  absorption  of  carton  dioxide  from  the  at¬ 
mosphere.  These  were  overcome  by  'using  freshly  boiled 
water  for  solutions  and  by  passing  nitrogen  over  any  vessel 
open  to  the  air  during  the  preparation  of  solutions. 

Results 


The  equation  for  the  reaction  studied  is  that  for 
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second-order— first-order  reversible  kinetics*: 

,t 


P  +  CN  <==*  PCN 
klj. 


(1) 


The  course  of  the  reaction  can  be  followed  conveniently 
using  light  of  wavelength  422  mp,  which  was  found  to  corre¬ 
spond  to  the  maximum  of  the  Soret  absorption  peak  of  the 
HRP-CN  complex  (this  value  differs  slightly  from  Keilin  and 
Hartree's  assignment  of  423  mP  for  the  same  maximum6 ) .  A 
photograph  of  a  typical  oscilloscope  trace  showing  a  record 
of  photomultiplier  output  voltage  vs.  time  is  reproduced  in 
Figure  14.  Values  of  these  quantities  were  read  from  an 
enlargement  of  the  picture  and  the  results  were  converted  to 
concentration  of  the  HRP-CN  complex  vs.  time  in  the  manner 
described  below. 

As  pointed  out  in  Chapter  2,  the  photomultiplier  out¬ 
put  voltage  is  proportional  to  the  transmittance  of  the  so¬ 
lution  in  the  stopped-flow  observation  chamber.  If  the 
reaction  is  slow  compared  to  the  mixing  time,  as  is  the  case 
with  the  reaction  of  cyanide  with  HRP,  one  can  assume  that 
the  point  on  the  trace  corresponding  to  t-0,  henceforth 
designated  TP,  represents  the  transmittance  of  pure  HRP  only 
since  effectively  no  reaction  has  occurred.  Similarly,  the 


The  symbols  i — »  are  used  for  reactions  occurring  at  measur¬ 
able  rates  and  - — N  for  reactions  in  which  equilibrium  is 


assumed  to  be  maintained 


6  o 


« 
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Figure  14.  Typical  oscilloscope  trace  showing  the 
P^oiT^oss  of  the  reaction  between  HRP  and  cyanide. 
Horizontal  scale:  500  mv/division;  vertical  scale: 
0.5  sec/ division . 
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voltage  output  when  the  reaction  reaches  equilibrium  corre¬ 
sponds  to  the  transmittance  of  the  reaction  mixture,  TE. 

One  can  measure  TP  simply  by  taking  a  portion  of  the  stock 
peroxidase  to  be  used  in  the  stopped-flow  experiment,  dilut¬ 
ing  by  a  factor  of  two,  and  measuring  the  transmittance  of 
the  solution  on  a  Beckman  DU  spectrophotometer.  Similarly, 
a  1:1  mixture  of  stock  HRP  and  ligand  solutions  can  be  meas¬ 
ured  on  a  DU  to  yield  TE.  In  both  cases,  a  small  correction 
must  be  made  to  take  into  account  the  fact  that  path  length 
in  the  DU  is  1.000  cm,  and  the  path  length  in  the  stopped- 
flow  apparatus  is  I.067  cm.  The  transmittance  at  any  time, 

T ( t ) ,  can  then  be  calculated  by  linear  interpolation.  If 
a  photograph  of  an  experimental  trace  is  projected  onto 
graph  paper  and  the  point  TP  is  set  equal  to  0,0  one  can 
obtain  numerical  readings  from  the  graph  paper,  R(t),  as  a 
function  of  time,  and  also  a  numerical  reading,  RE,  which 
corresponds  to  TE  for  the  reaction.  One  can  then  write  the 
interpolation  equation: 

T ( t )  =  TP  -  (TP“TE)  (2) 

The  values  of  T(t)  can  be  converted  to  corresponding  absorb¬ 
ance  values  by: 

a( t )  =  log[l.0/T(t)]  (3) 

We  now  have  values  for  total  absorbance  vs.  time  for  the 
course  of  the  reaction.  However,  the  total  absorbance  at 
422  mp  has  contributions  from  both  HRP  and  HRP-CN  and,  at 
any  time  during  the  course  of  the  reaction,  is  given  by: 


' 

1 


at  =  (£p[p]  +  £c[pcn]  )i  o  067  (4) 

where  £p  and  £  c  are  the  molar  absorptivities  of  HRP  and 
HRP-CN  at  422  myij  and  lo067  the  path  length  in  the  stopped- 
flow  observation  chamber <.  From  eg  1,  we  can  introduce  into 
eq  4  the  expression 

[P]  =  [P]o  -  [PCN]  (5) 

which  results  in 

a-t  =  1  Oo6?(£p[p]o  +  (£c  “£P)  [pcn]  )  (6) 

£c  was  found  to  be  constant  over  the  pH  range  used  in  this 
study j  however 9  £P  varied  considerably  at  pH  values  above 
9o 5  due  to  the  shift  in  the  Soret  peak  of  HRP  upon  formation 
of  its  hydroxide  complex «  Equation  6  can  be  simplified  with 
additional  auxiliary  absorbance  measurements 0  Another  por¬ 
tion  of  the  stock  HRP  solution  is  taken  and  added  to  an 
equal  volume  of  cyanide  solution  of  sufficient  concentration 
that  100  °/o  PCM  is  formed^  and  the  absorbance  of  the  solu¬ 
tion  is  measured  to  yield  a  term  designated  AC0  Knowing  AC 
and  AP  (the  absorbance  corresponding  to  TP)S  an  equation 
can  be  written  which  expresses  £p  in  terns  of  £t»  For  a 
path  length  of  loOO  cm? 


AP 

=  £P[p]o 

(7) 

and 

AC 

= 

(8) 

Equations 

7  and  8 

can  be  combined  to  gives 

€P 

_  A££s 
“  AC 

(9) 

Now ^  substituting  eq  7  and  9  into  eq  6  s 
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A-r  =  1.06?(AP  +  £t(l  -  ||)[PCN])  (10) 

Equation  10  can  be  rearranged  to  give  the  final  desired 

equation  which  gives  [PCN]  at  any  time: 

At  —  067AP  ( 1 1  ) 

1.0676.(1  - 

For  each  experiment  which  was  performed,  duplicate 
solutions  were  made  up  for  the  measurement  of  each  of  AP,  AC 
and  TEj  these  values  were  measured  and  the  average  of  the 
determinations  calculated.,  After  the  stopped-flow  experi¬ 
ment  was  performed,  values  of  RE  and  R(t)  vs.  time  were  read 
from  the  photographs  of  the  oscilloscope  traces  with  the  aid 
of  a  photographic  enlarger.,  Data  consisting  of  the  values 
of  AP,  AC,  TE,  RE  and  R(t)  and  the  corresponding  times, 
along  with  a  value  of  £c=  9*39  x  104  M“1ctTx  as  determined 
by  Keilin  and  Hartree6,  were  used  with  a  computer  program 
consisting  of  eq  2,  3>  an<3  11  to  yield  values  of  [PCN[]  vs. 
time.  Once  [PCN]  is  known  as  a  function  of  time  for  the 
course  of  the  reaction,  one  can  proceed  with  the  analysis  of 
the  data  from  the  photographs  of  the  oscilloscope  traces . 

The  differential  rate  equation  corresponding  to  eq  1  is: 

=  ki[P][CN]  -  k~ i [PCN ]  (12) 

The  apparent  rate  constants  were  obtained  by  direct  use  of 
eq  12,  but  with  d[PCN]/dt  approximated  byA[PCN]/&t;  i.e., 
by  taking  a  chord  over  a  short  interval  as  equal  to  the  tan¬ 
gent  at  the  center  of  that  intervale  Using  this 


'JA 
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approximation,  along  with  the  average  concentration  of  PCN 
complex  within  the  interval,  [PCN] eq  12  becomes: 

=  *([P]o  -  [PCN]aJ([CN]o  -  [PCN],,)  -  k'lCPCN^v; 

(13) 

In  order  that  the  chords  would  be  good  approximations  of 
tangents,  the  following  technique  was  used:  ^  [PCN]/^t  was 
determined  over  short  intervals  of  time  in  the  early  stages 
of  the  reaction  when  the  rate  was  changing  rapidly  with 
time,  and  then  at  increasingly  larger  intervals  as  the  reac¬ 
tion  approached  equilibrium.  Values  of  k'x  and  k-i  could 
then  be  obtained  from  eq  13  when  experimental  data  consist- 
ingA[PCN]/At  and  the  corresponding  [PCN]^  was  used  with  a 
nonlinear  least- squares  computer  program.  A  copy  of  this 
program  is  found  in  Appendix  I.  Fifteen  to  twenty  points 
from  each  experimental  trace  were  used  to  help  ensure  a  val¬ 
id  analysis.  Generally  three  oscilloscope  traces  were  ana¬ 
lyzed  for  each  experiment  and  average  values  of  k’i  and  k?-i 
were  determined. 

Although  eq  12  can  be  integrated  when  [p]o  ^  [CN]0» 
the  resultant  equation  is  complex49.  For  the  case  where 
[P]0  =  [CN]0,  the  integrated  form  of  eq  12  is61: 

[FCKK[CN]o3  -  [PCM][PCH])  =k.t  m) 
([CN]oS  -  [PCT7]a  )  [CN]oa([TCT] -[PCN]) 

A  plot  of  the  left  side  of  this  equation  vs.  time  should 
give  a  straight  line  of  slope  equal  to  ki  and  an  intercept 
of  zero.  Unfortunately,  eq  14  suffers  from  two  disadvan¬ 
tages:  it  does  not  yield  k-i  and  it  can  be  used  only  if 


. 
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[P]o  equals  [CN]0o  A  plot  of  the  left  side  of  eq  14  vs, 
time  for  an  experiment  at  pH  6 o5  where  [P]o  was  equal  to 
[CN]0  is  shown  in  Figure  15°  The  resulting  value  for  kJ  i 

agreed  within  5  °/o  of  the  value  obtained  using  the  computer 
analysis  of  eq  13  o 

The  experimental  rate  constants  obtained  at  each  pH 
studied  are  presented  in  Table  IIIo  Found  in  Table  IV  are 
values  of  the  dissociation  constants  for  the  HRP-CN  complex 
obtained  in  two  different  ways  s  direct  spectrophotometric 
determination  and  calculation  from  s  kLj,/k\o  The  spectro 
photometric  determination  consisted  of  determining  jjPCN]  for 
an  equilibrium  mixture  of  HRP  and  cyanide  and  using  eq  15 > 
below s  for  calculation  of  the  equilibrium  constant 0 


Plots  of  the  logs  of  k*i,  and  vs,  pH  are  shown  in 

Figures  16$  17  and  18$  respectively o 

Most  of  the  work  was  done  using  crude  HRP 5  runs  at 
pH  4o2s  5°5  and  80 5  were  repeated  using  pure  HRPo  The  re- 
suits  obtained  from  the  pure  HRP  agreed  within  10  °/o  of  the 
corresponding  results  using  crude  HRP0  As  will  be  discussed 
further  in  Chapter  5>  the  possible  presence  of  isozymes1 °> 1 1 
in  the  commercial  HRP  used  does  not  affect  significantly  the 
observed  kinetic  results  of  this  study  or  their  interpreta¬ 


tion  0 


■ 
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Figure  15.  Test  of  reversible  second-order— first-order 
kinetics.  The  left  side  of  eq  14  is  plotted  vs.  time 
using  data  for  the  reaction  of  HRP  with  cyanide  at  pH  6»5» 
The  slope  is  equal  to  k’i . 
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Table  III 

Rate  Data  for  the  Binding  of  Cyanide  by  Peroxidase  at  25  °C. 

Experimental  Rate  Constants 
With  Standard  Deviations 


pH  ki  (M”1sec“'1)  k'-i  (sec”1) 


4 

o  21 

(6 

.l±0o3) 

X 

104 

do 

7*0.2) 

X 

10“ 

4 

.46 

(7 

0*00) 

X 

104 

(1. 

7*0.2) 

X 

10“ 

4 

.  9  6 

(8 

. 8*0 0  5 ) 

X 

104 

(1. 

9*0.2) 

X 

10“ 

5 

.50 

(1 

0  2*0  el) 

X 

106 

(1. 

9*0.2) 

X 

10“ 

5 

.94 

(1 

0  0*0  el) 

X 

106 

(1. 

8*0 . 3 ) 

X 

10“ 

6 

.50 

(9 

o7±0o5) 

X 

104 

(2. 

0*0.2) 

X 

10" 

7 

.  05 

(9 

0  8*0 . 5 ) 

X 

104 

(2. 

8*0 . 3 ) 

X 

10" 

7 

.47 

(9 

•  0*0 . 5 ) 

X 

104 

(2. 

9*0.3) 

X 

10" 

8 

oil 

(7 

.2*0.4) 

X 

104 

(2. 

9*0.3) 

X 

10" 

8 

°55 

(6 

.6*0.4) 

X 

104 

(2. 

9*0.2) 

X 

10“ 

8 

.98 

(4 

.6*0.3) 

X 

104 

(3- 

2*0.2) 

X 

10" 

9 

09 

(2 

.4*0.2) 

X 

104 

(2. 

9*0.2) 

X 

10" 

9 

o84 

(1 

.0*0.1) 

X 

104 

(2. 

3*0.1) 

X 

10" 

10 

olO 

(5 

.9*1.1 ) 

X 

103 

(2. 

6*0.3) 

X 

10" 

10 

03 

(3 

.4*0.3 ) 

X 

103 

(1. 

9*0.2) 

X 

10" 

10 

.48 

(2 

.2*0.2) 

X 

103 

(2. 

0*0.2) 

X 

10" 

10 

03 

(1 

0  0*0 . 1 ) 

X 

103 

(1. 

3*0.1 ) 

X 

10" 

11 

.05 

(2 

.7*0.3) 

X 

10s 

(8. 

9*1.2) 

X 

10" 

11 

Ol 

(1 

.0*0.1 ) 

X 

102 

(7° 

3*1 . 0 ) 

X 

10" 

' 
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Table  IV 


Equilibrium  Data  for  the  Binding  of  Cyanide 
by  Peroxidase  at  25 °C. 


pH 


Kd  = 

-  [POT] 


From  Kinetic  Data 


Spectrophotometric 

Determination* 


4.21 

(2. 9 ±0.4) 

X 

10~6 

(2. 6*0.2) 

X 

10"6 

4.46 

(20*00) 

X 

10"6 

2.1 

X 

10"  6 

4.96 

(2. 1*0. 3) 

X 

10”6 

2.0 

X 

10“6 

5*50 

(1.7*0. 3) 

X 

10-6 

(1.5*0. 2) 

X 

10-6 

5.94 

(1.7*0. 3) 

X 

10~6 

1 . 6 

X 

10~6 

6.50 

(2. 1*0. 2) 

X 

10“  6 

(1.9*0. 2) 

X 

10-6 

7  •  05 

(2. 9*0. 3) 

X 

10-6 

(2. 4*0. 2) 

X 

10“6 

7.47 

(3. 2*0. 4) 

X 

10"6 

2.7 

X 

10-6 

8.11 

(4.1*00) 

X 

10“  6 

(30*0.1 ) 

X 

0 

1 

0) 

8.55 

(40*0.4) 

X 

10“6 

(4. 0*0. 2) 

X 

1  q—  6 

8.98 

(7*  0*0.5 ) 

X 

10“6 

(6.3*00) 

X 

10~6 

90? 

(1. 2*0.1) 

X 

10-5 

(1.2*0. 1) 

X 

10-5 

9.84 

(2. 3*0.1) 

X 

10~6 

2.2 

X 

10“  6 

10.10 

(4. 5*0.7) 

X 

10-6 

(4.2*0. 7) 

X 

10"  5 

1003 

(5. 7*0. 6) 

X 

10“6 

5.4 

X 

Hr 6 

10.48 

(9. 1*1. 4) 

X 

10~s 

( 8  0  8*0.4) 

X 

10"  6 

10.73 

(1.3*0. 2) 

X 

10-4 

(10*0.1) 

X 

10”4 

11.05 

(30*0.6) 

X 

10~4 

(3. 1*0. 4) 

X 

10“4 

11.31 

(7. 1*1. 2) 

X 

10-4 

6 . 6 

X 

10“4 

limits 

are  average 

deviations 

for  multiple 

determina- 

tions,  no  limits  given  for  single  determinations. 
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Figure  16.  Plot  of  log  ki  vs,  pH  for  the  binding  of 
cyanide  by  HRP .  The  solid  line  is  calculated  using  the 
best-fit  parameters  of  mechanism  V. 


*• 
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Figure  17 •  Plot  of  log  k-i  vs.  pH  for  the  dissociation  of 
HRP-cyanide  complex.  The  solid  line  is  calculated  using 
the  best-fit  parameters  of  mechanism  V. 


log  K 
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pH 


Figure  18.  Plot  of  log  Kd  vs.  pH  for  the  HRP-CN  complex. 
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Discussion 

The  pK  for  the  dissociation  of  hydrocyanic  acid  is 
9° 21  at  25 °C.  and  zero  fa63.  An  extended  Debye-Huckel  equa¬ 
tion63  was  used  to  calculate  the  activity  coefficients  at 
p  equal  to  0.11  and  yielded  a  corrected  pK  of  9*0*  Since 
the  study  was  made  over  a  pH  range  encompassing  both  sides 
of  this  pK,  both  CN"  and  HCN  must  be  considered  as  possible 
binding  ligands.  The  first  mechanism  considered  is  the 

binding  of  cyanide  to  HRP  with  no  heme-linked  ionizable 
* 

groups  : 

CN"  +  P  PL 

Kl  1  ,  ^ 
HCN  +  P  PL 

The  representation  of  mechanism  I  contains  several  simplify¬ 
ing  features:  hydrogen  ions  have  been  omitted,  the  charges 
on  the  proton  species  are  not  shown  since  they  are  unknown, 
and  the  HRP-ligand  complex  has  been  designated  PL  without 
specifying  whether  the  ligand  is  bound  as  HCN  or  CN"  since 
this  is  irrelevant  to  the  analysis  of  the  ki  data.  The  form 
of  the  complex  and  the  nature  of  the  reverse  reaction  will 

heme-linked  ionizable  group  is  defined  as  any  ionizable 
group  in  HRP  which  influences  the  binding  of  a  ligand  at  the 
sixth  coordination  position  of  the  heme  ferric  iron.  This 
group  may  be  part  of  either  the  prosthetic  group  or  protein 
moiety. 
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be  considered  for  the  analysis  of  the  k’-i  data. 

For  the  analysis  of  the  apparent  forward  rate  constant 
data,  it  is  necessary  to  have  an  equation  which  relates  the 
specific  rate  constants  of  mechanism  I  to  k’i  of  eq  1.  From 


mechanism  I,  neglecting  the  reverse  rate  term: 

=  ki[p][CN~]  +  ka[P][HCN]  (16) 

dt 

Multiplying  the  right  side  of  eq  16  by 
( [HCN]+[CN~ ] ) / ( [HCN]+[CN~ ] )  yields : 

d[PL]  =  kx[P][CN~]([HGN]+[CN~])  k8 [P ][HCN] ( [HCN>[CN~ ] )  (  } 

^t  [hcn]+[cn~]  [HCN]+[CN“] 

Introducing  the  notation  [CN]  =  [CN“]  +  [HCN]  and  rearrang¬ 


ing,  eq  17  becomes: 

<3[PL]  =  ld[P][CN]  t  Ka[P][CM] 

dt  [cm-]  i  [hcm]  [hcn]  [cn~ 

[CN“]  [CN~]  [HCN]  [HCN 


Simplification  of  eq  18  by  use  of  the  ligand  dissociation 
constant  KL=  ( [H+][CN" ])/[HCN]  leads  to  the  expression: 


a[PL] 

dt 


ki[p][CN] 

77HT 


K,. 


+ 


ka[P][CN] 

,  Ku 
1  +  T — Tn 


(19) 


Comparing  eq  19  with  eq  12  one  obtains  the  desired  result, 


which  is: 


k'i 


ki 


1  + 


[H+] 

Kl 


JSs. 


1  + 


KL 

[H+] 


(20) 


Equation  20  can  be  rearranged  to  the  form  of  a  linear 


equation : 


■ 
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(21) 


If  mechanism  I  were  valid,  a  plot  of  k*i(l  +  Kl/[h+])  vs. 


l/[H+]  would  be  linear,  with  a  slope  equal  to  kiKL  and  an 
intercept  equal  to  ka 0  As  can  be  seen  in  Figure  19s  this 
plot  is  not  linear,  so  mechanism  I  is  not  valido 

The  next  simplest  mechanism  involves  the  participation 
of  one  heme-linked  group; 


(II) 


Total  charges  on  the  peroxidase  species  and  the  equilibria 
involving  the  ligand  species  have  been  omitted,  although 
the  effect  of  the  dissociation  of  HCN  on  the  concentration 
of  ligand  has  been  taken  into  account  in  the  mathematical 
treatment  of  the  mechanism o  A  derivation  analogous  to  that 
used  for  eq  20  yields  for  mechanism  II: 


(22) 


Equation  22  can  be  used  with  a  nonlinear  least- squares  pro¬ 
gram  to  yield  the  best-fit  values  of  specific  rate  constants 
and  equilibrium  constant  for  the  heme-linked  group.  Because 
the  k"i  values  have  a  range  of  103,  the  nonlinear  least- 


IJ 
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Figure  19  •  Test  of  mechanism  for  HCN  and  CN“  binding  to 
peroxidase  with  no  heme-linked  acid  groups  (mechanism  i). 
If  this  mechanism  were  valid  the  plot  of  k'i  ( 1  +  Kl/[h+]) 
vs.  would  be  linear. 
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squares  program  was  modified  for  these  and  all  subsequent 
analyses  to  minimize  the  sum  of  the  squares  of  the  relative 
residuals  rather  than  the  sum  of  the  squares  of  the  absolute 
residuals *  Details  of  this  modification  are  found  in  Appen¬ 
dix  Io  If  k3  and  k5  are  set  equal  to  zero,  eq  22  represents 
a  mechanism  where  HCN  is  postulated  to  be  the  only  effec¬ 
tive  attacking  ligand,  and  if  k4  and  k6  are  set  equal  to 
zero  it  represents  CIST  as  the  only  effective  attacking  lig¬ 
and  <>  Although  all  combinations  were  tried  using  nonlinear 
least-squares  analysis,  it  was  not  possible  to  fit  the 
data  within  experimental  error*  Thus  mechanism  II  is  elim¬ 
inated* 

The  next  mechanisms  to  consider  involve  two  heme-linked 
ionizable  groups*  A  large  number  of  possible  combinations 
of  attacking  ligands  and  bound  ligands  can  be  conceived  for 
such  a  mechanism*  Only  two  of  those  possibilities  will  be 
briefly  discussed  -  the  only  two  which  came  close  to  fitting 
the  experimental  rate  data  using  nonlinear  least-squares 
analysis  * 
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Ka 


HP  + 
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P  + 


POH 


HCN 


HCN 


^  HPHCN 


k- 


Ka 


■  i  o 


k“io 


*  PHCN 


(IV) 


From  these  mechanisms,  one  can  write  the  following  equations: 

(23) 


A(1 

Kl 


— iix- _ _  +  — — - - _ _ 

+  JSfc  +  i+I£_J  +  _£a_ 


[H"0  '  [H+]a 


Ka  [H+] 


^  +I|^) 


ifi. 


JSUQl 


K; 


1  +  CH+] 

[H+]  [Ht]s  +  Ka  +  [H+] 


•j_  +  Ka  +  KaK3 


(24) 


t  6 

k~i 


JB. 


1  + 


Kac 

[H+] 


1  + 


is n 

Kac 


(25) 


k- !  =  ___ (2  6) 

1  +  ^ 

[H+]  Kac 

Equations  23  and  25  correspond  to  mechanism  III,  and  eq  24 
and  26  correspond  to  mechanism  IV 0  By  proper  choice  of  the 
values  for  the  adjustable  parameters  in  both  eq  23  and  24 , 
a  satisfactory  fit  to  the  k’x  data  is  obtained*  The  best-fit 
parameters  and  their  standard  deviations  obtained  from  non¬ 
linear  least-squares  analysis  are  shown  in  Table  V»  The 
specific  rate  constants  k8  and  kg  have  effective  values  of 


zero,  which  means  that  although  they  might  have  finite 
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Table  V 

Rate  and  Equilibrium  Constants 
Obtained  from  Analysis  of  Equations  23  and  24 

Mechanism  III,  Eq.  23 

k7  =  (1. 0*0.1)  x  lO10  M“1sec~1 
0  —  k8 ^  20  M“1sec”1 
K3  =  (1. 1±0.1)  x  10“4  M 
Ka  =  (1.5*0. 4)  x  10"11  M 


Mechanism  IV,  Eq.  24 

0  <  k9<  1  x  104  M’1sec“1 
kio  =  (9*6*0.?)  x  104  M_1sec“1 

t 

Ka  =  (1. 1*0.1)  x  10“4  M 
K3  =  (1.5*0. 4)  x  10~11  M 
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values,  their  corresponding  terms  in  eq  23  and  24  do  not 
affect  significantly  the  predicted  values  of  k'l.  They  are 
included,  however,  because  of  the  discussion  below  of  the 
principle  of  detailed  balancing .  Upper  limits  on  their 
values,  shown  in  Table  V,  are  estimated  by  assuming  that  the 
corresponding  terms  in  eq  23  and  24  contribute  a  maximum  of 
five  per  cent  to  the  predicted  values  of  k’i « 

The  predicted  value  of  the  pK  of  the  most  basic  heme- 
linked  group,  10.8,  appears  to  correspond  to  the  pK  for  the 
formation  of  the  HRP-hydroxide  complex.,  This  pK  has  been 
measured  at  the  values  ranging  from  10.6  to  11.319j34» 
Therefore  P  and  POH  have  been  used  in  mechanisms  III  and  IV 
to  represent  HRP  with  the  most  acid  heme-linked  groups  ion¬ 
ized,  and  with  water  and  hydroxide,  respectively,  in  the 
sixth  coordination  position  of  the  heme  ferric  iron. 

One  flaw  of  mechanisms  III  and  IV  is  that  they  allow 
no  simple  explanation  for  the  behavior  of  k-i  at  high  pH. 

For  the  pH  region  4o2-9°4,  nonlinear  least-squares  analyses 
of  eq  25  and  26  give  the  best-fit  parameters  and  their 
standard  deviations t  k»7  =  k~9  =  0. 17*0. 02  sec“l; 
k-e  s  k-io  =  0o30±0o02  sec"1;  K3c  =  (lo9±0c8)  x  10“7  M. 

Although  both  mechanisms  III  and  IV  can  account  for 
the  pH  dependence  of  the  k’i  data,  they  are  not  valid  be¬ 
cause  the  principle  of  detailed  balancing  is  violated.  By 
applying  this  principle  to  mechanism  III,  it  can  be  shown 


that ; 


■ 


k.7  6, K5  c 

k-  7K3 
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(27) 


k8  = 


Using  reverse  rate  constant  data  from  above  and  other  data 
from  Table  V  in  the  right  side  of  eq  27  >  one  obtains  a 
detailed  balance  value  of  k8  =  (3°0±1.6)  x  107  M^sec”1, 
while  -the  direct  kinetic  analysis  results  in  a  value  of 
ks  ^  20  M”1sec~1.  Similar  analysis  applied  to  mechanism  IV 
results  in  ks  =  (3*1*1 «6)  x  10 7  M“1sec““1  compared  to 
k9  ^  1  x  104  M”1sec”10  In  addition*  in  the  case  of  mech¬ 
anism  IV s  evidence  from  inorganic  chemistry  indicates  that 
cyanide  complexes  of  transition  metals  do  not  exist  with 
HCN  as  the  bound  ligand64'65.  If  mechanism  IV  is  revised 
so  the  ligand  is  bound  as  the  anion,  it  is  not  possible  to 
fit  the  k?-i  data  with  the  equation  derived  from  this  ver¬ 
sion.,  Thus  the  mechanisms  III  and  IV  must  be  eliminated. 

It  should  be  mentioned  here  that  the  protein  ioniza¬ 
tion  schemes  represented  in  mechanisms  III  and  IV  are  also 
simplifications.  That  ionization  scheme  can  be  more  proper- 
ly  written  as: 


HP  OH 


^  /  Kib 

POH 


(28) 


which  shows  that  there  are  really  four  ionization  constants 
and  four  species  involved,  rather  than  two  ionization  con¬ 
stants  and  three  species,  as  implied  in  mechanisms  III  and 
IV o  Following  the  treatment  of  Dixon  and  Webb66 


» 


Ks  and  K3 , 


' 
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the  experimentally  observable  constants,  can  be  considered 
"molecular"  ionization  constants  and  Kj,cl,  K2a.,  Kib  and  K2b 
the  true  group  ionization  constants.  If  the  two  ionizing 
groups  are  sufficiently  far  enough  apart  that  they  do  not 
influence  each  other,  then  Km  =  Kib  and  Ks  a.  =  K2b.  It  can 
be  shown66  that  the  relationship  between  the  molecular  ion¬ 
ization  constants  and  the  group  ionization  constants  iss 


(29) 


Ka  =  Kjlq.  +  K2 


(30) 


and 


If  Kicl  »  K2a,  then  eq  29  reduces  to  K2  =  K1Q_  and  if  Klb»K2b, 
then  eq  30  reduces  to  K3  =  K2b»  Thus,  although  the  kinetic 
study  yields  only  molecular  ionization  constants,  if  these 
constants  differ  by  more  than  a  factor  of  100,  they  can  be 
taken  to  represent  group  ionization  constants  quite  accu¬ 
rately. 

Similarly,  since  the  rigorous  treatment  shows  that 
there 'are  four  possible  species  present  rather  than  the 
three  indicated  by  the  simplified  treatment,  the  meaning  of 
the  kinetic  assignments  must  be  examined.  The  binding  of  a 
ligand  species  to  HP  and  POH  are  unambiguous,  but  what  was 
represented  as  binding  of  a  ligand  species  to  P,  shown  in 
mechanisms  III  and  IV,  could  actually  be  binding  to  either 
or  both  of  the  species  P  or  HPOH  represented  in  eq  28.  This 
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ambiguity  cannot  be  resolved  by  use  of  the  kinetic  method* 
The  rate  expressions  for  the  binding  to  the  species  P  and 
HP OH  from  eq  28  can  be  written: 

Rate  of  binding  to  P  =  ka[p][L]  (31) 

Rate  of  binding  to  HPOH  =  kb[HPOH][L]  (32) 

where  [l]  represents  the  ligand  concentration.  From  the 
equations  for  the  group  ionization  constants  it  can  be 
shown  that : 

[p]  =  (33) 

[HPOH]  K3  a. 

Substituting  [HPOH^J  from  eq  33  into  eq  32 : 

Rate  of  binding  to  HPOH  =  kb(|^)[P][L]  (34) 

Equation  3^  is  of  the  same  form  as  eq  31j  tut  with  ko_ 
replaced  by  kb(Kaa,/K1(x )  •  Thus,  it  is  impossible  to  tell 
whether  the  rate  constants  k8  and  ki  0  correspond  to  ko., 
kb>  or  a  linear  combination  of  the  two. 

We  finally  consider  mechanism  V  which  involves  three 
heme-linked  ionizable  groups: 


s  1 
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k- 
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?= 
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(V) 


K5 ,  Ke,  K4c  and  K5C  are  molecular  ionization 
,  The  corresponding  rate  equations  are: 


ki  i  +  + 


j Si  5K4il5 


fH+l  ’  [H+]3 


( 1  .  K4  K4K 6  K4KsK6  Ch  .] 1 

IF]  [H+]»  [H+]'  m  Kl  > 


j.  ki 4K4  ,  ki fiK4Ks 
ki  a  +  Tvffi  + 


H 


/  „  K*  K4K5  KaKbKb  \  ,  Kt 

1  [h^]  T&y  [h+] 
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i  = 


lr  Th+I  -i.  jSr  1  3K4C  i_  rr  ,  fe~Ji J^acK ,5tfc  4.  KriL 

+  k'laLH  J  +  TH^]  +  k-l*K*<-  +  TH+r  +  H+T 


1  +  T& 


KicK^i 


[ip]  [H+] 


(36) 


Four  simplifications  of  mechanism  V  and  of  the  correspond¬ 
ing  eq  35  and  36  fit  the  experimental  rate  data  satisfacto¬ 
rily  and  do  not  violate  the  principle  of  detailed  balanc¬ 
ing.  These  simplified  versions  of  mechanism  V  all  contain 
two  pairs  of  forward  and  reverse  specific  rate  constant 
terms.  ,  All  other  terms  do  not  contribute.  The  forward 
rate  constants  for  these  four  simplifications  are: 

(A)  kn,  ki  3 ;  (B)  ki4,  kls;  (C)  ki3>  ki4;  (D)  kn,  ki6. 

The  various  kinetic  and  equilibrium  parameters,  obtained 
from  nonlinear  least-squares  analysis,  are  listed  in  Tables 
VI  and  VII.  The  same  equilibrium  constants  for  the  heme- 
linked  ionizable  groups  and  identical  predictions,  to  four 
significant  figures,  of  the  ki  and  k’-i  data  are  obtained 
from  the  four  mechanisms.  The  values  of  k*i  and  k-i  computed 
from  the  best-fit  parameters  of  the  four  mechanisms  are 
shown  by  the  solid  lines  in  Figs.  16  and  17 »  respectively. 

The  equivalence  of  the  four  mechanisms  in  their  fit  to 
the  experimental  rate  data  is  a  result  of  symmetry  proper¬ 
ties  inherent  in  the  overall  mechanism  V.  Thus  the  rate 
of  binding  of  CN"  by  HSP,  given  by  ki i[HaP][CN“],  can  be 
rearranged  using  the  relations  K*  =  [h+][HP]/[hsp]  and 
Kl  =  [H+][CN"]/[HCN]  to  give  (kiiKL/K*)[HP][HCN].  Similarly 


’ 
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Table  VI 

Rate  and  Equilibrium  Constants 
of  Simplified  Versions  A  and  B  of  Mechanism  V 


Version 


Constant 

A 

B 

kii 

(7.441.3) 

X 

109 

M“1sec“1 

* 

kis 

* 

* 

kia 

(3. 8*0.4) 

X 

107 

M“1sec”1 

* 

ki4 

* 

(lol±0c 

>2) 

X 

10B 

M^sec”1 

ki  5 

* 

# 

ki6 

(9-3*1 . 

,0) 

X 

104 

M”1sec”i 

k-  1 1 

0.17=fcOo03 

•“  i 

sec 

k-i3 

* 

k-i3 

0»30i0o03 

sec” 1 

k-u 

* 

(9.4*1. 

>6) 

X 

106 

M^sec”1 

k_  i  5 

* 

* 

k- 1 6 

# 

(1.4*0. 

.1) 

X 

101 

0  M"1sec”1 

k4 

( 7  o  2  £2  o  0 ) 

X 

10” 

5  M 

(7.2*2. 

.6) 

X 

10” 

6  M 

Ks 

(4.3*0°6 ) 

X 

10“ 

7  M 

(4.3*4. 

>8) 

X 

10” 

7  M 

Ks 

(lo4±0»4) 

X 

10“ 

ii  M 

(1.4*0, 

>4) 

X 

10“ 

11  M 

(1  <.8*0.9 ) 

X 

10“ 

7  M 

(  1  o  8*0  < 

>9) 

X 

10“ 

7  M 

Ksc 

(2. ,1*0. 8) 

X 

10” 

ii  M 

( 2  o  1  *0 . 

»8) 

X 

10” 

11  M 

*  signifies  a  rate  constant  effectively  equal  to  zero 


' 
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Table  VII 

Rate  and  Equilibrium  Constants 
of  Simplified  Versions  C  and  D  of  Mechanism  V 

Version 


Constant 

C 

D 

fcn 

(7o4±1o3) 

X 

109 

M”1 sec” 1 

kia 

% 

* 

ki  3 

( 3 » 8±0 

•  4) 

X 

107 

M”1sec“’1 

ki4 

(1 ol^O 

.2) 

X 

10B 

M"1sec“1 

* 

ki  s 

# 

* 

kie 

* 

(9. 3*1-0) 

X 

104 

M”1sec”1 

k-n 

* 

0cl7i0.03 

sec"1 

k-i  3 

* 

* 

k- 13 

o .30^0 

o03 

sec”1 

k-u 

(9 0^*1 

06) 

X 

10s 

M“1sec“1 

* 

k-  i  s 

* 

k-°  1 6 

* 

( 1 o4±0 0 1 ) 

X 

101 

0  M”1 sec”1 

k4 

(?c2±2 

*  6 ) 

X 

10” 

5  M 

(7o2*2o0) 

X 

10" 

6  M 

k6 

(403*0 

0  6 ) 

X 

10“ 

7  M 

( 4  0  3  *4  •  8 ) 

X 

10" 

7  M 

k6 

(1.4±0 

.4) 

X 

10" 

11  M 

(Io4±0o4) 

X 

10” 

u  m 

K*c 

(lo8±0 

.9) 

X 

10" 

7  M 

(lo8±0.9) 

X 

10” 

7  M 

Ksc 

(2ol±0 

.8) 

X 

10“ 

11  M 

( 2  0 1&0 .8) 

X 

10" 

11  M 

*  signifies  a  rate  constant  effectively  equal  to  zero 


' 
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ki3[HP][CN"]  can  be  converted  to  (kj.3Ku/Ke)[p][HCN]. 

The  equivalent  fit  of  mechanisms  VA  and  VB  to  the  experi¬ 
mental  kinetic  data  implies  that  ki*  =  knKL,/K4,  and 
ki6  =  ki3KL/Ke.  These  relations  are  obeyed  by  the  data  in 
Table  VI®  In  a  similar  manner,  it  can  be  shown  from  the 
kLi  data  that  k-14  =  k-.11/K4c  and  k-i6  =  k-i3/K6Ci.  It 
follows  that  the  four  mechanisms,  VA  -  VD,  are  kinetically 
indistinguishable®  It  also  follows  that  appropriate  linear 
combinations  of  the  four  mechanisms  will  provide  equally 
satisfactory  fits  to  the  experimental  rate  data® 

As  an  example  of  the  principle  of  detailed  balancing, 
it  can  be  shown  that : 


ki3 


k,ijLk-.i.a£*&. 

k- 1 1K4 


(37) 


From  values  from  Table  VI  for  the  factors  in  the  right  side 
of  eq  37 »  ki3  =  (3®3=tl  -9 )  x  107  M”1sec“1  compared  to 
ki3  =  (3°8=±0<.4)  x  107  M”1sec“1  obtained  from  direct  kinetic 
analysis  o 

In  addition  to  the  four  simplified  versions  of  mechan- 
ism  V  discussed  above,  one  other  mechanism  appears  to  fit 
the  experimental  rate  data®  It  is  similar  to  mechanism  VB 
but  involves  the  formation  of  a  complex  with  undissociated 
HCN  as  the  bound  form  of  the  ligand®  Aside  from  evidence 
from  inorganic  chemistry,  it  can  be  eliminated,  as  was  the 
case  with  mechanism  IV,  by  applying  the  principle  of 
detailed  balancing® 


_ 
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The  four  simplifications  of  eq  35  which  fit  the  experi¬ 
mental  data  can  in  turn  be  simplified  in  certain  pH  regions , 
depending  on  the  relative  magnitudes  of  the  various  para¬ 
meters  o  These  simplified  equations  can  be  rearranged  into 
linear  form  to  provide  an  alternative  means  of  calculating 
some  of  the  parameters,  as  a  check  on  the  nonlinear  least- 
squares  analysis.  For  example,  in  the  pH  region  4 .21-8 .11, 
the  expression  for  ligand  binding  according  to  mechanism  VA 
can  be  reduced  to: 


k'di  +I2l4(i  +  ^ 


+ 


Ks 


)  =  hi  3  + 


Kn[H+] 


kl  k4  [h+]'  ”  k4 


(38* 


A  plot  of  the  left  side  of  eq  38  [H+]  for 

K4  =  2  x  10” 6  M  and  K6  =  4.3  x  10“ 7  M  is  shown  in  Fig.  20. 

Linear  least-squares  analysis  of  the  plot  yields  the  values 
kn/K4  =  (1.0±0o2)  x  1014  sec”1  and 

ki3  =  (3o6=t0.6 )  x  107  M“1sec“1o  Similarly,  for  the  pH 
region  8o55”ll°31>  eq  35  reduces  to: 


kid  +  l5U)(i  +  =  ki 

KL  [H+]  [H+] 


(39) 


Using  the  values  K5  =  4. 3  x  10~7  M  and  K6  =  1.4  x  10”11  M, 
calculation  of  the  left  side  of  eq  39  yields  a  mean  value 
of  ki 3  =  (3.9*0. 3)  x  107  M“ 1  sec” 1 . 

The  predicted  pK  value  for  the  most  acid  heme-linked 
group,  4.1,  is  consistent  with  results  of  other  studies38’41 
The  pK  value  of  6.4  does  not  agree  with  values  of  5 • 0  and 
7.0  reported  elsewhere38119  but  is  in  agreement  with  a  re¬ 
analysis  of  the  HRP-fluoride  kinetics67.  The  value  of 
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Figure  20.  Plot  of  the  left  side  of  eq  38  vs.  [H+] 
used  to  test  the  validity  of  mechanism  VA  for  the  bind¬ 
ing  of  cyanide  by  HRP  over  the  pH  range  4.2  to  8.1. 
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this  intermediate  pK  obtained  from  the  HRP-cyanide  kinetic 
data  is  a  sensitive  parameter.  For  example,  it  is  changed 
from  a  value  of  6.4  to  7»2  by  arbitrarily  changing  the  val¬ 
ue  of  pKL  from  9 « 0  to  9 • 1 «  On  the  other  hand  the  pK  values 
of  4.1  and  10.8  are  only  changed  by  0.1  for  the  same  change 
in  pK  Lo 

The  values  of  the  equilibrium  constants  of  the  heme- 
linked  groups  are  shifted  to  smaller  values  as  a  result  of 
the  binding  of  the  ligand.  This  shift  is  in  the  direction 
expected  if  the  bound  ligand  is  negatively  charged68. 

Large  values  of  kii  and  k_i6  are  obtained  from  mechan¬ 
ism  V,  raising  the  question  of  whether  these  values  might 
exceed  the  diffusion-controlled  limit,  which  for  a  second- 
order  specific  rate  constant  for  an  enzyme-ligand  reaction 
is  of  the  order  of  1010  M~1sec“1  .  Within  the  error  lim¬ 
its  of  the  various  parameters  required  for  computing  the 
diffusion-controlled  limit,  ku  and  k«i6  do  not  exceed  this 
value.  Therefore,  mechanisms  VA,  VB  and  VD  cannot  be  elim¬ 
inated  on  this  basis. 

In  summary,  the  cyanide  kinetic  results  indicate  the 
presence  of  three  heme-linked  groups  with  pK  values  of  4.1, 
6.4  and  10.8  on  horseradish  peroxidase.  They  also  indicate 
that  the  anion  form  of  the  ligand  exists  in  the  peroxidase- 
ligand  complex,  but  it  is  impossible  to  distinguish  which, 
if  either,  fqrm  of  the  ligand  dominates  the  binding  reac¬ 
tion.  The  same  ambiguity  exists  as  to  the  reactivities  of 
the  various  peroxidase  species. 
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CHAPTER  4.  THE  EFFECT  OF  LIGAND  BINDING  AND  ACID  SPLITTING 

ON  THE  OPTICAL  ROTATORY  DISPERSION  OF  HRP 

Introduction 

A  beam  of  plane-polarized  light  may  be  considered  to 
be  made  up  of  two  circularly  polarized  rays,  one  whose  vec¬ 
tor  rotates  clockwise  and  one  counterclockwise  as  the  beam 
advances.  A  medium  containing  asymmetric  molecules  trans¬ 
mits  the  two  components  with  unequal  velocity,  i.  e, ,  it 
has  different  refractive  indices  for  the  two  circularly  pol¬ 
arized  rays.  If  neither  ray  is  absorbed,  the  two  rays  have 
a  resultant  which  is  a  plane-polarized  ray,  but  whose  vi¬ 
bration  is  in  a  plane  which  has  been  rotated  with  respect 
to  that  of  the  incident  ray«  !The  plane  of  polarization  of 
the  resultant  moves  through  a  steadily  increasing  angle  as 
the  beam  passes  through  the  medium <>  The  amount  of  rotation 
is  measured  as  the  angle  between  the  incident  and  emergent 
beams,  and  is  normally  expressed  as  specific  rotation,  [a]A9 
which  is  given  bys 

Mr 

where  is  the  angular  rotation  of  the  sample  in  degrees 
using  light  of  wavelength  A  9  1  is  the  path  length  in  deci¬ 
meters,  and  w  is  the  number  of  grams  of  solute  in  one  ml  of 
solution.  Since  the  indices  of  refraction  are  dependent  on 
wavelength,  different  wavelengths  of  light  are  rotated 
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different  amounts,  giving  rise  to  the  phenomenon  known  as 
optical  rotatory  dispersion** 

If  the  medium  also  absorbs  the  two  circularly  polar¬ 
ized  components  unequally,  the  emergent  light  beam  is  el- 
liptically  polarized  and  the  medium  exhibits  circular  di- 
chroism.  The  combined  phenomena  of  optical  rotation  and 
circular  dichroism  occurring  in  the  wavelength  region  cor¬ 
responding  to  an  optically  active  absorption  peak  is  known 
as  a  Cotton  effect*  The  ORD  curve  of  a  Cotton  effect  is  a 
symmetrical  curve  shaped  like  a  derivative  of  an  absorption 
peak  and  is  centered  at  a  wavelength  corresponding  to  the 
maximum  of  the  absorption  peak*  The  maximum  of  a  Cotton  ef 
feet  curve  is  known  as  the  peak,  and  the  minimum  as  the 
trough.  When  the  peak  occurs  at  a  wavelength  longer  than 
that  of  the  maximum  of  the  absorption  curve,  the  Cotton  ef¬ 
fect  is  termed  "positive”,  and  when  it  is  at  a  shorter  wave 
length,  the  Cotton  effect  is  "negative" * 

Cotton  effects  have  been  classified  into  two  groups. 
Using  the  terminology  suggested  by  Blout70,  "intrinsic  Cot¬ 
ton  effects"  reflect  the  optical  rotatory  power  of  the  pro¬ 
tein  proper,  and  "extrinsic  Cotton  effects"  are  generated 

■^Abbreviations  used  in  this  chapters  ORD,  optical  rotatory 
dispersion;  CD,  circular  dichroism;  HRP,  HRP-CN,  HRP-F , 
HRP-OH,  HRP-Na ,  ferric  horseradish  peroxidase  and  its 
cyanide,  fluoride,  hydroxide,  and  azide  complexes,  respec¬ 
tively;  PGA,  poly-<x-L-glutamic  acid. 


-78- 


by  a  local  chromophoric  site  on  the  protein  or  by  a  chromo- 
phoric  molecule  interacting  with  an  asymmetric  site  on  the 
protein*  A  good  review  of  the  use  of  extrir^sic  Cotton 
effects  for  the  study  of  interactions  between  proteins  and 
biologically  important  molecules  such  as  coenzymes,  sub- 
strates,  inhibitors,  and  prosthetic  groups  containing  metals 
has  been  written  by  Ulmer  and  Vallee71. 

In  I96I  Simmons,  et  al.73  proposed  that  the  trough  at 
233  mji  in  the  ORD  spectra  of  many  proteins  was  conformation 
dependent  and  that  its  depth  provided  a  measure  of  the 
°t-helical  content  of  the  protein*  Although  there  is  some 
debate  as  to  whether  one  can  calculate  accurately  the  ab¬ 
solute  amount  of  ^-helical  structure  using  Simmons’  meth¬ 
od73’74,  it  is  certainly  a  sensitive  tool  for  measuring 
changes  in  protein  conformation* 

Recently  this  method  has  been  used  to  study  the  effect 
of  ligand  binding  on  the  conformation  of  heme  proteins*  It 
has  been  found  that  while  myoglobin  and  catalase  show  no 
change  in  helical  content  upon  the  binding  of  ligands 73 ' 7 6 > 76 , 
and  while  hemoglobin  shows  no  change  upon  the  binding  of 
most  ligands77,  hemoglobin  seems  to  undergo  an  8*5  °/o  de¬ 
crease  in  helicity  upon  binding  of  oxygen78*  A  similar 
study  carried  out  on  HRP  would  be  expected  to  show  if  it 
behaves  like  most  cases  involving  heme  proteins,  or  whether 
a  gross  conformational  change  takes  place  upon  the  binding 
of  ligands*  Any  conformational  changes  which  might  occur 
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when  the  prosthetic  group  is  split  from  the  protein  could 
also  be  observed  by  an  ORD  study. 

Ulmer  and  Vallee79  have  shown  from  a  study  of  the  ORD 
of  HRP ,  catalase,  and  their  cyanide  complexes  that  the  ex¬ 
trinsic  Cotton  effects  apparent  in  the  Soret  region  reflect 
the  orientation  of  the  heme  with  respect  to  the  protein. 

Their  work  has  been  extended  to  include  the  fluoride,  azide, 
and  hydroxide  complexes  of  HRP  and  HRP  solutions  in  which 
the  heme  is  partially  or  completely  removed  from  the  protein. 

Experimental 

Analytical  reagent  grade  HRP  (RZ  =  2»9>  lot  6485206) 
was  obtained  from  the  Boehringer-Mannheim  Corp. ,  N.  Y. ,  and 
was  dialyzed  against  water  and  centrifuged  before  use. 

Hemin  of  White  Label  grade  was  obtained  from  Eastman  Organic 
Chemicals.  Inorganic  chemicals  were  reagent  grade  and  were 
used  without  further  purification.  HRP  concentrations  were 
determined  from  absorbance  measurements  at  403  and  497  mp 
using  a  Beckman  DU  spectrophotometer  and  literature  values 
for  the  molar  absorptivities 6 .  Typical  concentrations  of 
HRP  used  were  2.5  x  10“ 6  M  for  ORD  measurements  in  the  UV 
region  and  1.5  x  10" 5  M  for  measurements  in  the  visible 
region.  All  pH  measurements  were  made  with  a  Beckman  ex¬ 
panded  scale  pH  meter.  ORD  spectra  were  obtained  at  room 
temperature  using  a  Jasco  Model  ORD/UV-5  recording  spectre- 
meter  set  at  a  scan  rate  of  1.7  mp/min. 
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Absorption  spectra  of  solutions  were  obtained  using  a 
Perkin  Elmer  model  202  recording  UV-visible  spectrophotom¬ 
eter.,  For  all  solutions  used  in  this  study,  the  maximum 
absorbance  was  less  than  two  in  the  wavelength  regions  of 
interest o  This  procedure  was  followed  to  help  ensure  that 
the  ORD  spectra  obtained  were  not  caused  by  rotatory  ar¬ 
tifacts  produced  by  regions  of  high  absorbance80.  The  val¬ 
ues  of  specific  rotation  were  calculated  using  the  molar 
concentrations  determined  spectrophotometrically  and  a 
value  of  40,200  for  the  molecular  weight  of  HRP12. 

The  HRP-F  complex  was  studied  at  pH  4*3  in  a  solution 
which  was  10~2  M  with  respect  to  KF;  HRP-CN  at  pH  6.0  with 
10"4  M  KCN ;  HRP-OH  at  pH  11. 3  obtained  by  addition  of  KOH; 
and  HRP-N3  at  pH  4.3  with  0.5  M  NaN3 .  From  the  dissociation 
constants  of  the  complexes1 6 » 41 > 6 7 ,  it  was  estimated  that 
the  HRP-CN,  HRP-F,  and  HRP-N3  complexes  represented  greater 
than  95  °/o  of  the  total  HRP  in  solution,  and  that  HRP-OH 
was  about  75  °/o  formed.  The  limiting  factor  for  HRP-OH 
formation  is  the  instability  of  the  protein  at  high  pH30. 

Techniques  for  splitting  the  heme  group  from  the  pro¬ 
tein  moiety  of  HRP  have  been  described  in  a  series  of  pa¬ 
pers  by  Maehly30”32.  For  this  study,  HRP  was  split  by 
making  the  solution  acidic  and  allowing  it  to  react  at  room 
temperature  for  several  hours.  The  splitting  reaction,  fol¬ 
lowed  by  taking  absorption  spectra  of  the  solution30,  was 
not  complete  in  2 6  hours  at  pH  3°0»  but  appeared  to  be 


. 
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complete  after  12  hours  at  pH  2o 75 »  No  precipitate  of  heme 
or  denatured  protein  was  observed  during  the  course  of  the 
splitting  experiments. 

A  solution  of  hemin  was  prepared  in  an  ethanol-water 
mixture  containing  44.5  weight  percent  ethanol,  a  medium  in 
which  hemin  is  soluble  and  stable.  The  final  solution  was 
1.5  x  10" 5  M  in  hemin,  and  h?  x  10"4  M  in  NaOH. 

Results 

Figure  21  shows  the  ORD  curves  for  HRP  and  split  HRP 
from  215  rnp  to  350  mp,  plotted  as  specific  rotation  vs. 
wavelength.  Within  experimental  error,  the  ORD  curves  of 
the  hydroxide,  fluoride,  and  cyanide  complexes  of  HRP  are 
the  same  as  that  of  pure  HRP.  The  ORD  of  HRP-N3  could  not 
be  obtained  in  this  region  due  to  the  high  absorbance  of 
azide.  The  curves  in  Figure  21  show  a  minimum  at  233  nip 
which  is  characteristic  of  the  conformation-dependent 
trough  in  the  negative  intrinsic  Cotton  effect  centered  at 
225  rnp.  For  a  quantitative  measure  of  the  amplitude  of 
this  trough,  the  reduced  mean  residue  rotation,  [R'jaaas 

73 

has  been  used,  as  suggested  by  Simmons,  et  al.  The  values 
of  [R* ]s33  for  HRP,  its  cyanide,  fluoride,  and  hydroxide 
complexes  and  for  split  HRP  were  calculated  from  the  formula ; 

1 

=  3M55-]aa£-T  (2) 

100 (na33  +  2) 

where  the  refractive  index  for  water,  n33s$  was  taken  as 
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Figure  21.  ORD  curves  of  HRP  ( - )  and  split  HRP  (-  -  -)  in  the  UV  region 
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l*39sl>  and  the  mean  residue  molecular  weight,  Mr,  as  11580. 
It  is  reasonable  to  assume  that  the  rotation  observed  at 
233  caused  by  protein  only.  Since  there  is  18  weight 

percent  carbohydrate  material  present  in  HRP10>X3,  a  cor¬ 
rected  specific  rotation,  [oc:]333,  was  used  in  eq  2,  based 
on  an  effective  protein  molecular  weight  of  40,200  x  0.82. 
Table  VIII  summarizes  the  results  for  HRP,  its  complexes 
and  split  HRP,  where  the  tabulated  value  of  [R5]a33  repre¬ 
sents  the  mean  and  the  error  limit  the  mean  deviation  of 

the  determinations*,  Multiple  determinations  were  carried 

✓ 

out  for  two  reasons  s  to  ascertain  the  reproducibility  of 
the  ORD  spectrum  for  a  given  sample,  and  to  check  on  the 
variation  of  [R'jaas  with  change  in  concentration  of  the 
protein.  Although  the  variation  in  HRP  concentration  was 
greater  than  a  factor  of  two,  small  mean  deviations  of 
[R*3a33  were  obtained.  In  addition,  the  [R5]a33  data  show¬ 
ed  no  trend  with  concentration;  only  a  random  scatter  about 
the  mean  was  observed.  The  lack  of  concentration  dependence 
of  [r 1 J333  is  indicative  that  the  Cotton  effect  observed  is 
genuine80  0 

Figure  22  shows  the  extrinsic  Cotton  effects  associ¬ 
ated  with  the  heme  Soret  bands  of  HRP  and  its  fluoride, 
cyanide,  and  hydroxide  complexes.  Similarly,  Figure  23 
shows  the  ORD  curves  for  HRP-Na,  completely  split  HRP,  and 
HRP  which  was  partially  split  by  reacting  it  at  pH  3.0  for 
20  hours.  The  ORD  curve  of  the  1.5  x  10" 6  M  hemin  solution 
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Table  VIII 


Reduced  mean  residue  rotation  at  233  for 
horseradish  peroxidase,  its  hydroxide,  cyanide, 
fluoride  and  azide  complexes,  and  for  split  HRPo 


Species 

-[R']s»a  x  10"3 

Number  of 
Determinations 

HRP 

6 o  0^0  »2 

7 

HRP- OH 

6  - 0±0o 1 

3 

HRP-CN 

5«9±0.1 

3 

HRP-F 

5<,9=t0ol 

3 

split  HRP 

5  # 1 *0  o 1 

2 
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Figure  22.  ORD  curves  in  the  visible  region  for  HRP 
and  its  fluoride,  cyanide  and  hydroxide  complexes. 
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Figure  23.  ORD  curves  in  the  visible  region  for  the 
azide  complex  of  HRP  and  for  HRP  which  has  had  the  heme 
group  partially  and  totally  split  from  the  protein. 
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was  examined  in  the  region  365-^50  mji#  The  optical  rotation 
was  constant  in  that  wavelength  region,  with  no  evidence  of 
a  Cotton  effect  at  the  Soret  band  of  the  hemin<> 

Discussion 

Although  [R'Jsaa  for  the  HRP-ligand  complexes  appears 
smaller  than  for  pure  HRP,  the  difference  is  within  exper¬ 
imental  error*  This  same  result  was  observed  in  the  case 
of  myoglobin  and  its  complexes73*  It  can  therefore  be  con¬ 
cluded  that  the  binding  of  ligands  to  HRP  has  less  than 
2  °/o  effect  on  the  amount  of  ^-helical  conformation  of  the 
protein. 

Splitting  the  heme  from  the  protein  portion  of  HRP 
results  in  a  15  °/o  decrease  in  [R'^Jaas®  This  result  agrees 
with  the  difference  found  between  the  values  of  [R*]aaa 
for  apomyoglobin  and  metmyoglobin7 6 > 83 ,  and  supports  the 
suggestion  that  the  heme  binding  stabilizes  the  molecule  as 
a  whole20.  Inspection  of  Figure  21  shows  that  the  ORD 
curve  of  split  HRP  is  lower  than  that  of  HRP  in  the  region 
250-320  mja,  a  result  which  was  not  observed  in  similar 
studies  on  myoglobin  and  apo-myoglobin7 6 1 82 ,  but  comparison 
of  the  ORD  curves  of  catalase  and  acid-denatured  catalase 
shows  the  same  phenomenon76.  This  would  indicate  that  for 
HRP,  there  may  be  other  Cotton  effects  from  the  heme  group 
which  affect  the  ORD  in  this  region,  and  although  no  such 
fine  structure  was  observed  in  the  present  study,  results 
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from  a  recent  study  of  the  CD  of  HRP  showed  the  presence  of 
bands  at  282,  335,  and  370  m p83. 

The  value  of  [R* ]aaa  =  -6000  for  HRP  indicates  that  it 
has  less  helical  content  than  myoglobin,  catalase,  or  hemo¬ 
globin,  all  of  which  have  values  of  [R*]aa3  of  about 
-90007 6  * 76  ' 7a>  83  and  that  HRP  has  a  greater  helical  content 
than  cytochrome  c,  which  has  a  value  of  [R*]aaa  of  about 
-440084,  Use  of  [R’jaaa  to  calculate  the  helical  content 
of  proteins  seems  to  result  in  an  underestimate  in  many 
cases80,  which  may  be  a  result  of  the  fact  that  fully  coil¬ 
ed  PGA  has  been  used  as  the  standard  for  100  °/0  helicity. 
Although  there  seems  to  be  agreement  that  for  random  PGA, 
[R*]aaa  =  ”2000,  there  is  a  wide  disparity  in  the  reported 
values  for  fully  coiled  PGA73 ' 74 > 76 •  Further,  there  seems 
to  be  an  optical  anomaly  associated  with  what  had  been  sup¬ 
posed  to  be  the  fully  coiled  form8 6 > 86  <>  Thus,  it  appears 
that  a  better  estimate  of  the  helical  content  of  HRP  can 
be  made  using  the  value  of  £R*]aaa  for  myoglobin  as  a  sec¬ 
ondary  standard,  since  X-ray  crystallographic  studies  have 
shown  that  it  contains  77  °/o  <X-helix87  1  88 «  Using  [R']a3a 

73  ,76)78,82,8© 

equal  to  -2000  for  0  °/o  helix  and  -9200  for 

77  °/o  helix,  and  assuming  that  the  heme  Cotton  effects 
make  a  negligible  contribution  to  [R1 ]a33S  HRP  is  calcu¬ 
lated  to  contain  43  °/o  CL-helical  structure,  which  is  re¬ 
duced  to  33  °/o  for  split  HRP,  This  estimate  of  the  helical 
content  of  HRP  agrees  with  a  value  of  40  0 / o  helical  content 
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found  in  a  study  of  the  CD  of  HRP83 . 

The  ORD  curves  for  HRP  and  HRP-CN  in  the  visible  region 
agree  with  those  found  previously79.  The  curves  for  HRP-F, 
HRP-OH,  and  HRP-N3  are  also  of  the  same  sign  and  order  of 
magnitude  as  HRP  and  HRP-CN,  thus  leading  to  the  conclusion 
that  the  orientation  of  the  heme  with  respect  to  the  pro¬ 
tein  is  the  same  for  HRP  and  its  ligand  complexes.  While 
the  ORD  curve  of  HRP  in  the  visible  region  is  very  similar 
to  those  of  hemoglobin  and  myoglobin89,  it  differs  signif¬ 
icantly  from  those  of  catalase79  and  ferricytochrome  c84>90>91 
Further,  the  ORD  curves  of  HRP-CN,  ferricytochrome  c-cyanide91 
and  f errimyoglobin-cyanide73  are  all  similar  to  each  other, 
but  as  a  group  are  quite  different  from  that  of  catalase- 
cyanide79.  It  is  known  that  the  heme  iron  of  cytochrome  c 
has  histidine  and  methionine  bound  in  the  fifth  and  sixth 
coordination  positions  and  that  hemoglobin  and  myoglobin 
have  histidine  and  water  in  the  fifth  and  sixth  positions92. 

A  comparison  of  absorption  spectra,  magnetic  moments,  and 
reduction  potentials  of  the  heme  proteins  led  Brill  and 
Williams  to  conclude  that  HRP  "is  an  amino- carboxy late,  or 
less  probably  a  water-amino  complex,  and  that  catalase  is  a 
dicarboxylate  or  less  probably  a  carboxy late-water  com¬ 
plex"23.  However,  Brill  and  Sandberg  later  concluded  that 
for  HRP,  histidine  is  in  the  fifth  position27.  The  results 
of  the  present  study,  when  compared  with  the  ORD  curves  of 
the  other  heme  proteins  as  mentioned  above,  indicate  that 
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HRP  is  probably  similar  to  myoglobin  and  hemoglobin,  with 
histidine  in  the  fifth  position  and  water  in  the  sixth 
position. 

The  results  on  HRP-OH  disagree  with  an  earlier  find¬ 
ing  that  upon  the  formation  of  alkaline  HRP ,  the  ordered 
structure  of  the  protein  is  destroyed93.  The  results  of  the 
present  studies,  both  at  233  m p  and  in  the  Soret  region, 
indicate  little  difference  between  the  structures  of  HRP 
and  HRP-OH,  Osbahr  and  Eichhorn  reached  their  conclusion 
from  ORD  measurements  of  an  extrinsic  Cotton  effect  associ¬ 
ated  with  the  absorption  peak  at  640  mp  in  HRP.  This  Cotton 
effect  disappeared  upon  raising  the  pH  from  7*2  to  11.0, 
a  result  which  can  alternatively  be  interpreted  as  being  due 
to  the  disappearance  of  the  absorption  peak  at  640  mp  when 
HRP-OH  is  formed 6 • 

Inspection  of  Figure  23  shows  that  as  the  heme  splits 
from  the  protein,  the  positive  Cotton  effect  in  the  visible 
region  broadens  and  diminishes  in  magnitude,  and  then  dis¬ 
appears  altogether  when  the  heme  has  been  completely  split 
from  the  protein.  This  result,  along  with  the  observed 
lack  of  any  Cotton  effect  for  the  hemin  solution,  confirms 
the  contention  that  the  extrinsic  Cotton  effect  arises  from 
the  interaction  of  the  heme  with  an  asymmetric  site  on  the 
protein79 . 

In  a  recent  study  on  the  inactivation  of  HRP  when 
treated  with  azide16,  Brill  and  Weinryb  concluded  that 
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azide  was  attacking  a  methionine  residue  at  the  active  site. 
Evidence  from  the  present  study  suggests  instead  that  the 
observed  inactivation  was  due  to  the  splitting  of  the  heme 
from  the  protein.  Their  experiment  in  which  HRP  was  in¬ 
activated  by  reacting  it  with  azide  at  pH  5.2  and  42°  for 
five  hours  was  repeated  in  this  study.  After  removal  of 
the  azide  by  dialysis,  the  Soret  absorption  spectrum  and  the 
ORD  curves  of  the  solution  at  both  the  UV  and  Soret  regions 
were  found  to  be  the  same  as  those  of  partially  split  HRP. 
Although  they  ran  a  control  in  the  absence  of  azide  at  pH 
5.2  and  found  no  inactivation,  it  is  known  that  splitting 
occurs  faster  at  a  given  pH  in  the  presence  of  ligands  such 
as  fluoride,  cyanide,  and  azide33.  Finally,  they  reported 
recovery  of  some  of  the  activity  of  the  inactivated  samples 
after  dialysis.  This  can  be  explained  on  the  basis  that  the 
splitting  is  reversible,  with  the  recombined  HRP  exhibiting 
normal  peroxidatic  activity33. 
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CHAPTER  5.  THE  EFFECTS  OF  TEMPERATURE  AND  IONIC  STRENGTH 
ON  THE  FORMATION  OF  THE  HYDROXIDE  COMPLEX  OF  HRP 

Introduction 

It  is  widely  accepted  that  HRP,  like  hemoglobin  and 
myoglobin,  contains  water  in  the  sixth  coordination  position 
of  the  heme  ferric  iron6^19'94,  although  the  suggestion  has 
been  made  that  the  heme  of  HRP  lies  in  a  crevice  with  both 
the  fifth  and  sixth  positions  occupied  by  protein  residues36 « 
When  a  solution  of  HRP  is  made  alkaline,  there  are  accom¬ 
panying  changes  in  the  spectrum  of  HRP6,  and  its  magnetic 
moment  changes  from  5*^5  to  2.66  Bohr  magnetons34,  indi¬ 
cating  the  formation  of  a  new  species,  HRP-OH*.  The  pK 
for  hydroxide  formation  has  been  measured  to  be  10*9  by 
spectrophotometry  and  11. 3  t>y  use  of  a  magnetic  balance34, 
and  10.6  by  a  study  of  the  reduction  potentials  of  HRP19, 
but  these  studies  were  not  carried  out  at  constant  ionic 
strength.  A  study  of  the  kinetics  of  cyanide  binding  to 

^Abbreviations  used  in  this  chapters  P,  refers  specifically 
to  a  form  of  HRP  in  which  it  is  assumed  that  water  is  in 
the  sixth  coordination  position;  HRP-OH  or  POH,  hydroxide 
form  of  HRP;  P,  POH,  H+,  equilibrium  values  of  P,  POH,  and 
H+;  Kc,  apparent  equilibrium  constant  defined  in  terms  of 
concentrations ,  and  hence  not  a  true  thermodynamic 
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HRP  at  an  ionic  strength  of  0.11  yielded  a  value  of  10.8  for 
this  pK  (Chapter  3)»  This  chapter  contains  a  report  on  a 
spectrophotometric  study  of  the  formation  of  HRP-OH  at  25»0° 
for  seven  values  of  ionic  strength  over  the  range  0.02  to 
0.20  and  also  at  18.0°  and  35»0°  with  p  =  0.11. 

Experimental 

Two  grades  of  HRP  were  obtained  from  Boehringer- 
Mannheim  Corp.,  N.  Y.:  ’’analytical  reagent  grade”  (RZ  =  2.9 
lot  6485206)  and  ’’purified  lyophilized  powder”  (RZ  =  0.8). 

All  HRP  was  exhaustively  dialyzed  against  conductivity  water 
and  centrifuged  before  use.  Analytical  reagent  grade  KN03 
(Mallinckrodt )  and  KOH  (B.D.H. )  were  used  without  further 
purification. 

Some  of  the  ’’purified”  HRP  was  treated  by  passage 
through  a  1.1  x  50  cm  column  of  C-50  CM  Sephadex  cation  ex¬ 
change  gel.  The  HRP  was  eluted  using  pH  6.2  phosphate  buf¬ 
fer,  with  the  ionic  strength  increasing  from  0.1  to  0.2 
during  the  course  of  the  elution.  A  single  passage  through 
the  column  resulted  in  the  recovery  of  a  fraction  of  RZ  = 

2.7  which  represented  65  °/o  of  the  initial  ’’purified”  HRP 
as  determined  by  absorbance  measurements  at  497  mju.  This 
purification  procedure  is  described  in  detail  in  Appendix  II. 
Absorbance  measurements  were  made  with  a  Beckman  DU  spec¬ 
trophotometer  which  had  the  cell  compartment  thermostated 
to  ±  0.1°.  A  Beckman  Expandomatic  pH  meter  was  used  for 


■ 
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all  pH  measurements. 

For  a  typical  experiment,  an  appropriate  amount  of 
1  M  KN03  and  enough  stock  HRP  to  make  a  final  HRP  concen¬ 
tration  of  about  1.7  x  10“ 6  M  were  pipetted  into  a  10  ml 
volumetric  flask,  which  was  then  filled  to  the  mark  with 
water.  The  10  ml  of  solution  was  transferred  to  a  30  nil 
beaker  which  was  covered  with  Parafilm  and  placed  in  a  con¬ 
stant  temperature  bath.  A  Beckman  39183  combination  pH 
electrode,  a  Teflon  syringe  needle,  and  a  platinum  stirring 
wire  were  inserted  through  slits  in  the  Parafilm.  The  pH 
was  measured  and  was  found  typically  to  have  a  value  of  5»5« 
A  three  ml  aliquot  of  the  solution  was  transferred  with  a 
syringe  to  a  cuvet  and  the  HRP  concentration  was  determined 
by  measuring  the  absorbance  of  the  solution  at  497  mp  where 
the  molar  absorptivity  coefficient  is6  1.00  x  104  M"1  cm"1. 
The  absorbance  of  the  solution  was  then  measured  at  4l6  m p, 
the  wavelength  corresponding  to  the  maximum  of  the  Soret 
peak  of  HRP- OH 5 o  A  second  HRP  solution  of  about  pH  5°5  was 
prepared  to  use  as  a  blank  in  the  spectrophotometer.  Its 
concentration  was  such  that  the  difference  in  absorbances 
between  the  blank  and  the  experimental  solution  was  about 
0.1  at  416  mp,  thus  keeping  the  absorbance  readings  during 
the  course  of  the  titration  in  the  range  0.1  -  0.6,  which 
includes  the  region  of  minimum  error  due  to  instrumental 
uncertainty69.  The  experimental  HRP  solution  was  returned 
to  the  beaker  and  the  pH  was  raised  by  the  addition  from 
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a  10  pi  Hamilton  syringe  of  a  few  pi  of  approximately  1  M 
KOH,  an  amount  small  enough  that  no  volume  correction  was 
necessary.  The  solution  was  stirred,  the  pH  was  remeasured, 
the  absorbance  of  an  aliquot  was  measured  at  4l6  mp  in  a 
capped  cuvet,  and  the  aliquot  was  returned  to  the  beaker. 
Nine  to  eleven  points  were  obtained  this  way,  with  the  ti¬ 
tration  being  carried  out  to  about  pH  11  **35°  Above  this  pH, 
the  heme  splits  from  the  protein  at  an  appreciable  rate30. 
After  each  increase  in  pH,  the  beaker  was  flushed  gently 
with  nitrogen  so  that  there  would  be  no  change  in  pH  or 
ionic  strength  due  to  absorption  of  carbon  dioxide  from  the 
atmosphere,  a  problem  encountered  particularly  above  pH  10. 
This  method  was  effective  as  proved  by  the  pH  stability  of 
the  solution  at  any  given  pH  over  a  period  of  15  minutes. 

While  it  is  true  that  the  addition  of  KOH  to  the  solu¬ 
tion  contributes  to  the  ionic  strength,  this  contribution 
is  negligible  for  most  cases.  At  pH  11.3,  the  KOH  added 
results  in  an  increase  in  ionic  strength  of  0.002.  This 
means  that  the  greatest  effect  due  to  the  added  base  would 
be  a  10  °/o  increase  in  the  ionic  strength  for  the  last 
point  taken  in  the  experiment  performed  at  p  =  0.02. 

Temperature- jump  kinetic  experiments  were  carried  out 
on  the  apparatus  which  was  described  in  Chapter  2.  A  tem¬ 
perature  jump  of  about  6°  was  applied  to  a  5° 18  x  10" 6  M 
solution  of  "analytical  reagent  grade"  HRP  at  25°  with 
pH  =  11.0  and  p  =  0.11.  The  photomultiplier  output  voltage 
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corresponding  to  a  change  in  absorbance  at  4l6  m \i  with  time 
was  displayed  on  an  oscilloscope *  Estimation  of  a  maximum 
relaxation  time  for  the  equilibration  of  the  HRP  -  HRP-OH 
mixture  was  made  by  examination  of  photographs  of  the 
oscilloscope  traces* 

If  it  is  assumed  that  a  water  molecule  is  in  the  sixth 
coordination  position  of  the  heme  iron  of  HRP j  its  ion¬ 
ization  can  be  represented  bys 

P  »  POH  +  H+  (1) 

Charges  on  the  protein  species  have  been  omitted  since  they 
are  not  known*  The  reversibility  of  the  process  represented 
in  eq  1  was  tested  in  the  following  experiment.  The  ab¬ 
sorbance  curve  of  the  Soret  peak  of  an  HRP  solution  at  pH 
5*5  was  taken,  the  pH  raised  to  11,  then  readjusted  to  pH  5°5 
and  the  spectrum  was  remeasured* 

Results 

The  reversibility  of  hydroxide  formation  was  demon¬ 
strated  by  the  invariance  of  the  spectrum  of  HRP  at  pH  5-5 
within  instrumental  error  (<0 *5  °/o )  before  and  after  adjust¬ 
ment  of  the  pH  to  11. 

The  total  absorbance  of  the  solution  at  any  pH  is  the 
sum  of  the  absorbances  due  to  the  P  and  POH  species,  which 
for  a  1  cm  cell  is  given  bys 

a  =  ep[p]  +  e  PO»  [poh]  (2) 

where  £p  and  £P0Hare  the  molar  absorptivities  of  P  and  POH 


. 
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at  4l6  myio  The  value  of  £  P0H  cannot  be  determined  directly 
because  it  is  not  possible  to  raise  the  pH  of  the  solution 
high  enough  to  form  100  °/o  POH  without  causing  irreversible 
changes  in  the  enzyme.  The  equilibrium  constant  corre¬ 
sponding  to  eq  1  is  given  by: 


Kc 


(3) 


Equations  2  and  3  can  be  combined  to  give  the  expression: 


A  A  =  [P]0A£-  £g!3*» 


(4) 


where  £A  =  A  -  A0,  where  A0  is  the  absorbance  of  the  solu¬ 
tion  when  all  HRP  is  in  the  P  form?  [P]o  is  the  initial 
concentration  of  HRP;  and  A£=  £P0 «  -  £P»  The  value  of  [H+] 
is  obtained  from  the  equations 

pH  =  ~log[H+]  (5) 

The  validity  of  eq  5  will  be  discussed  later. 

From  eq  4,  it  follows  that  a  plot  of  A A  vs«  [H+]AA 
should  yield  a  straight  line,  with  a  slope  of  -1/KC  and  an 
intercept  of  [p]0A£c  Such  a  plot  is  shown  in  Figure  24,  the 
data  being  from  an  experiment  using  ’’analytical  reagent 
grade”  HRP  at  25.0°  and  ju  -  0.14.  It  is  obvious  from  Fig¬ 
ure  24  that  the  experimental  data  do  not  fit  eq  4.  The 
deviation  from  eq  4  is  indicative  that  at  least  a  second 
process  is  occurring.  It  was  found  that  the  deviation  from 
the  expected  straight  line  was  smallest  when  ’’analytical 
reagent  grade”  HRP  was  used,  greater  for  Sephadex-treated 
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0  0.5  1.0  1.5 


AA[H+]  x  1011 

Figure  24.  Plot  of  AA  vs.  AA[H+]  for  a  spectrophotometric 
titration  of  pure  horseradish  peroxidase  at  25*0°  and 
p.  =  0.14.  The  lack  of  linearity  can  be  interpreted  to 
mean  that  more  than  one  ionization  is  taking  place. 
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"purified"  HRP  and  greatest  for  "purified”  HRP. 

From  the  temperature- jump  experiments  it  was  not  pos¬ 
sible  to  determine  exactly  the  relaxation  time  constant,  T* , 
for  the  reaction  being  observed  since  the  relaxation  appeared 
to  be  instantaneous  on  the  fastest  time  scale  which  could 
be  used.  The  time  scale  is  limited  by  the  instrument  re¬ 
sponse  time  determined  by  the  photomultiplier  load  resistor 
setting,  which  must  be  large  enough  to  obtain  sufficient 
sensitivity  to  observe  the  absorbance  change  accompanying 
the  chemical  relaxation.  However,  examination  of  photo¬ 
graphs  of  the  oscilloscope  traces  indicated  that  T/  must  be 
less  than  30  psec  using  the  solution  and  experimental  con¬ 
ditions  described  in  the  previous  section, 

v 

Discussion 

Among  the  possible  explanations  which  could  account  for 
the  observed  deviation  from  linearity  shown  in  Figure  24, 
one  of  the  simplest  might  be  that  a  heme-linked  ionizable 
group  with  a  pK  roughly  in  the  range  9“ 10  is  present  on 
HRP  and  that  its  ionization  gives  rise  to  another  spectro¬ 
scopically  distinguishable  species.  However,  this  possi¬ 
bility  is  ruled  out  since  the  amount  of  deviation  from 
linearity  of  plots  of  AA  vs.  AA£h+^J  varies  with  the  nature 
of  the  preparation.  A  more  plausible  explanation  is  that 
the  observed  effect  might  be  due  to  the  presence  of  more 
than  one  isozyme  of  HRP  in  the  preparations  used,  and  that 
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the  amount  of  the  minor  constituent  varies  with  the  prep¬ 
aration1  0  >  1 1  <, 

If  there  are  two  isozymes  present  which  have  the  same 
absorption  spectrum  for  the  P  species  but  which  have  differ¬ 
ent  equilibrium  contents  for  the  ionization  of  P  to  form 
POH,  one  can  derive  an  equation  analogous  to  eq  4  s 


a  A  = 


&£i[p]of 

77zp 

Kc  i 


A6»0]o(l-f) 


1  + 


~H+ 


] 


K 


ca 


(6) 


where  f  is  the  fraction  of  HRP  present  as  isozyme  1,  the 
subscripts  indicate  isozyme  1  or  isozyme  2,  and  the  other 
symbols  have  the  same  meaning  as  before 0  Equation  6  can 
be  used  with  a  nonlinear  least-squares  computer  program 
(Appendix  I)  to  yield  best-fit  values  of  A€if,  Kci,  A£a(l-f), 
and  Kca«  The  value  of  f  cannot  be  separated  from  that  of 

By  use  of  eq  6,  it  is  possible  to  obtain  a  good  fit  to 
the  experimental  titration  data  as  shown  in  Figure  25®  The 
circles  in  Figure  25  show  experimental  values  of  aa  as  a 
function  of  pH  while  the  solid  line  is  the  theoretical 
curve  obtained  from  use  in  eq  6  of  the  computer  best-fit 
values  of  KCi  =  3®?  x  1CT9  M,  A€j,f  =  5°3  x  10s  M“ 1  cm"1, 

Kca  =  9»1  x  10”la  M,  and  A£a(l~f)  =  3®9  x  104,  M"1  cm”1* 

Table  IX  summarizes  the  results  of  the  computer  anal¬ 
ysis  of  the  experimental  titration  data  using  eq  6o  These 
results  were  all  obtained  using  "analytical  reagent  grade” 
HRP  o  The  error  limits  shown  on  the  results  are  either 
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Figure  25.  Plot  of  AA  vs.  pH  for  a  spectrophotometric 
titration  of  pure  horseradish  peroxidase  at  25*0°  and 
ju  =  0.14.  The  solid  line  is  calculated  using  computer 
best-fit  parameters  in  eq  6. 


Summary  of  Results  of  Computer  Analysis  of 
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Duplicate  experiments  performed  for  these- .conditions . 
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standard  deviations  obtained  from  the  computer  analysis,  or 
in  the  case  of  duplicate  runs,  the  mean  deviation,  although 
for  the  duplicate  runs  the  two  types  of  errors  were  about 
equal*  The  values  of  A€if  andA£a(l“f)  show  no  variation 
with  change  of  temperature  or  ionic  strength;  the  mean  and 
mean  deviation  of  all  the  values  found  for  these  parameters 
ares  A6if  =  (5*4  ±  1*6)  x  10s  M"1  cm"1  and  A6a(l“f)  = 

(3*7  ±  0*1)  x  104  M”1  cm”1.  The  error  limits  quoted  for 
A£if  and  Kci  are  considerably  larger  than  those  for  A63(l-f) 
and  KCa,  which  is  a  consequence  of  the  relative  contribu¬ 
tions  of  the  two  terms  on  the  right  side  of  eq  6.  Although 
the  two  terms  are  typically  about  equal  at  pH  8.5?  the  sec- 
ond  term  is  ten  times  larger  than  the  first  at  pH  10. 5 »  and 
50  times  larger  at  pH  11  *3*  Thus  it  takes  a  relatively 
large  change  in  trial  values  of  A  6 if  or  Kci  to  influence 
the  overall  fit  of  eq  6  to  the  experimental  data. 

The  existence  of  several  isozymes  of  HRP  has  been  known 
for  some  time  7”9,  and  recently  the  properties  of  the  indi¬ 
vidual  isozymes  have  been  studied  in  detail10 >xlo  Shannon 9 
et  al.  reported  that  two  types  of  isozymes,  which  they  des¬ 
ignated  as  A~1  and  C,  had  different  spectral  properties10. 
They  showed  the  spectra  of  the  two  isozymes  at  pH  values 
ranging  from  4.4  to  10*5 •  In  neutral  or  acid  solution,  the 
Soret  absorption  bands  of  the  two  isozymes  are  almost  iden¬ 
tical.  The  spectrum  of  the  C  isozyme  at  pH  9.6  is  the  same 
as  that  near  pH  7°  However,  the  maximum  of  the  Soret  peak 
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of  the  A-l  isozyme  at  pH  9»6  shows  a  shift  to  higher  wave¬ 
length,  and  its  overall  absorption  spectrum  appears  to  be 
typical  of  that  of  partially  formed  HRP-OHo  At  pH  10 ,5 
both  isozymes  display  spectra  which  are  characteristic  of 
HRP-OH6. 

On  the  basis  of  these  findings,  it  appears  probable 
that  the  Boehringer-Mannheim  "analytical  reagent  grade"  HRP 
used  for  this  study  consists  of  a  minor  component  of  the 
A-l  type  of  isozyme,  and  a  major  component  of  the  C  type 
isozyme,  corresponding  to  the  designations  "isozyme  1"  and 
"isozyme  2"  used  in  eq  6.  Assuming  that  the  spectra  of  the 
P  and  POH  species  of  the  two  isozymes  are  the  same,  the 
relative  amounts  of  the  two  isozymes  present  can  be  esti¬ 
mated  from  the  mean  values  of  A£if  andA6a(l-f).  If  a£i 
and  A£a  are  equal,  then  (l-f)/f  =  3*7  x  104/5»4  x  102  and 
f  =  0.014,  thus  indicating  that  1.4  °/o  of  the  "analytical 
reagent  grade"  HRP  used  was  isozyme  of  the  A-l  type,  and  the 
rest  was  of  the  C  type.  This  calculation  also  allows  deter¬ 
mination  of  the  value  of  the  molar  absorptivity  coefficient 
of  HRP-OH  at  4l6  mp.  The  ratio  A403/A416  was  measured  for 
an  "analytical  reagent"  HRP  solution  at  pH  7  s  and  the  re¬ 
sult  0.668  was  obtained,  which  yielded  a  value  of  Cp4ie  « 
6.08  x  104  M“ 1  cm”1  from  Keilin  and  Hartree‘s  value6  of 
£p4oa  =  9  ol  x  104  M” 1  cm"1.  From  A^a  =  £PCH4i6  -  £p4is  = 
3.74  x  104  M"1  cm"1,  a  value  of  £PoH4ie  =  9.82  x  104  M_1  cm"1 
is  then  obtained,  which  differs  form  the  previously  reported 
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value6  of  £P0H4i6  =  8085  x  104  M"1  cm"1 «  However,  Keilin 
and  Hartree  determined  their  value  from  a  spectrum  of  HRP  at 
pH  11.4,  where  the  formation  of  HRP-OH  is  not  complete* 

Data  from  an  experiment  using  "purified  grade”  HRP 
indicates  that  there  is  about  12  °/o  A-l  type  isozyme  pres- 
ent  in  that  preparation*  Since  most  of  the  kinetic  study 
of  the  binding  of  cyanide  to  HRP  described  in  Chapter  3  was 
carried  out  using  "purified”  HRP,  the  question  is  raised  as 
to  whether  the  presence  of  two  isozymes  would  affect  the 
interpretation  of  the  results*  It  was  reported  that  exper¬ 
iments  repeated  at  three  pH  values  ranging  from  4 *2-8 *5 
using  "analytical  reagent"  HRP  gave,  within  experimental 
error,  the  same  observed  binding  rate  constants  as  those 
from  experiments  using  "purified"  HRP*  Therefore,  over  the 
pH  range  of  the  duplicate  experiments,  the  presence  of  two 
isozymes  did  not  affect  the  kinetic  results  since  the  rate 
constants  did  not  vary  with  the  relative  amounts  of  the  two 
isozymes  present*  At  higher  pH  values,  the  experimental 
data  from  the  cyanide  binding  study  was  fit  by  a  mechanism 
in  which  it  was  postulated  that  the  rate  of  binding  of  cya¬ 
nide  to  POH  was  negligible*  The  amount  of  POH  present  at 
any  pH  is  dependent  on  the  pK  for  its  formation,  which  the 
present  study  has  shown  to  be  different  for  the  two  isozymes* 
The  computer  best-fit  value  for  the  equilibrium  constant 
for  POH  formation  obtained  in  the  cyanide  study  was 
(1*4  ±  0*4)  x  1CT11  M  at  25°  and  p  =  0.11,  while  under  the 


' 
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same  conditions  and  using  the  same  definition  of  pH,  the 
present  study  yields  the  values  Kci  =  (lo8  ±  0.4)  x  10" 9  M 
and  Kcs  =  (1.00  ±  0.05)  x  10"11  M.  Thus,  the  equilibrium 
constant  for  hydroxide  formation  from  the  kinetic  study, 
which  according  to  the  present  interpretation  must  be  a 
weighted  mean  value  for  the  two  isozymes,  agrees  within  ex¬ 
perimental  error  with  Kcs  for  the  predominant  isozyme  2.  As 
a  further  check  of  possible  errors  in  the  previous  analysis, 
the  nonlinear  least- squares  analysis  of  the  forward  rate 
equation, given  by  equation  35  in  Chapter  3»  was  repeated  with 
terms  included  which  allowed  for  12  0 / o  isozyme  of  the  A-l 
type.  This  meant  that  the  term  containing  one  equilibrium 
constant  for  POH  formation  was  replaced  with  two  terms  con¬ 
taining  as  fixed  constants  the  values  of  KCi  and  Kcs  at 
p  =  0.11  found  in  the  present  study.  The  usual  iterative 
procedure  then  found  the  best-fit  values  for  the  other  rate 
and  equilibrium  constants.  Those  constants  were  all  found 
to  agree  within  5  °/0  of  the  values  found  in  the  previous 
study,  so  that  the  presence  of  two  isozymes  in  the  prep¬ 
aration  used  does  not  affect  the  previous  interpretation  of 
the  results  within  experimental  error. 

Although  a  precise  relaxation  time  was  not  determined 
for  the  reaction  involving  the  formation  of  POH,  it  is 
possible  to  establish  limits  for  the  rate  constants  in¬ 
volved.  Since  the  experiments  were  performed  using  "ana¬ 
lytical  reagent  grade"  HRP,  the  effect  of  the  small  amount 
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of  isozyme  A-l  will  be  ignored*  If  the  reaction  is  written 


as  an  associations 


POH  +  H+  P 

k;^ 

then  the  reciprocal  relaxation  time,  1/7"',  is  given  by63  s 


(7) 


=  ki([P0H]  +  [h+] )  +  k_i  (8) 

A  lower  limit  for  ki  can  be  obtained  from  eq  8  by  introduc¬ 
ing  the  relationship  k~i  =  kiKci,  and  the  values  of 
[POH]  =  2*6  x  lO"6  M,  [H+]  =  10"11  M,  and  T<  30  psec* 
Therefore,  30  x  104  sec”1  <  kj.(2*6  x  10"6  M) ,  or 
ki>l«3  x  lO10  M”1  sec"1  at  31°.  This  result  compares  with 
values  found  by  Eigen  and  co-workers  for  other  protonation 
reactions,  such  as  the  combination  of  a  hydrogen  ion  with 
a  hydroxide  ion,  which  has  a  rate  constant  of  1*3  x  1011  M”1 
sec"1  9S,  or  protonation  of  anions  such  as  formate,  benzoate, 
substituted  benzoates  and  phenolates,  all  of  which  have  rate 
constants  in  the  range  2  x  lO10  M"1  sec"1  to  6  x  lO10  Mml 


sec"1  at  25' 


9  6,97 


If  the  reaction  is  written  as  a  displacement  by  hydrox¬ 
ide  ions 

P  +  OH"  POH  (9) 

k-a 

introduction  of  Kw=  10" 14  M3  and  analogous  treatment  to 
that  above  allows  the  calculation  that  ka  >  1*6  x  107  M  1 
sec"1.  The  value  of  k3  appears  to  be  quite  large  in  com¬ 
parison  with  the  specific  rate  constants  for  fluoride  and 
cyanide  binding  to  HRP  with  the  two  most  acid  heme-linked 
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groups  ionized67,  and  although  the  process  represented  by 
eq  9  cannot  be  ruled  out  on  this  basis,  the  apparent  agree¬ 
ment  of  ki  with  specific  rate  constants  for  other  protonation 
reactions  is  evidence  that  the  proton  transfer  process  is 
the  one  being  observed  in  this  study* 

The  effect  of  ionic  strength  on  the  values  of  Kci  and 
Kcs  can  be  treated  on  the  basis  of  the  theory  of  the  sec¬ 
ondary  salt  effect.  For  a  reaction  of  the  types 

HAZ  *=>  H+  +  A(z“^  (10) 

where  Z  is  the  charge  on  the  undissociated  acid,  the  concen¬ 
tration  equilibrium  constant,  Kc,  is  related  to  the  thermo¬ 
dynamic  equilibrium  constant,  KT,  by  the  expressions 

log  Kc  =  log  KT  +  2(1  -  Z)D  (11) 

The  symbol  D  in  the  above  equation  stands  for  a  portion  of 
the  extended  Debye-Huckel  expression  due  to  Davies63,  which 
at  25°  and  with  water  as  the  solvent  is  given  bys 

D  =  0*509(7-^7-  "  0.20  (a)  (12) 

The  experimental  data  were  analyzed  using  eq  6  and  by  assum¬ 
ing  that  eq  5  is  valid.  However,  the  definition  of  pH  is 
an  operational  one.  The  term  pH  does  not  have  a  simple 
meaning;  it  is  not  a  true  measure  of  either  the  hydrogen 
ion  concentration  or  hydrogen  ion  activity,  but  probably  of 
a  quantity  between  the  two93.  If  it  were  a  measure  of 
hydrogen  ion  activity,  then  for  this  study,  eq  8  would  take 


the  forms 
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log  Kob*  =  log  Kt  +  (1-2Z)D  (13) 

where  K0bs  is  the  same  as  Ke  except  that  hydrogen  ion  concen¬ 
tration  is  replaced  by  hydrogen  ion  activity. 

Figures  2 6  and  27  are  plots  of  log  KCi  vs.  D  and  log  Kcs 
vs o  D,  respectively.  There  is  too  much  scatter  in  the  data 
in  Figure  26  to  try  to  interpret  quantitatively  the  results 
for  isozyme  1.  However,  inspection  of  Figure  26  shows  that 
the  plot  of  log  Kci  vs.  the  ionic- strength-dependent  term 
goes  through  a  minimum;  this  behavior  is  opposite  to  that 
found  in  corresponding  studies  on  the  formation  of  ferri- 
myoglobin  hydroxide"  and  f errihemoglobin  hydroxide100. 

The  data  shown  in  Figure  27,  for  isozyme  2,  fits  a  straight 
line  with  a  slope  of  -2.9±0.4  as  determined  by  linear  least- 
squares  analysis.  This  slope  may  be  used  to  predict  the 
effective  charge  on  the  P  species,  but  the  accuracy  of  this 
prediction  depends  on  the  interpretation  of  pH  and  hence  on 
whether  eq  11  or  eq  13  is  more  valid.  If  eq  11  were  obeyed, 
the  slope  of  Figure  27  would  be  expected  to  be  -2  for  Z  =  +2 
and  -4  for  Z  =  +3,  while  if  eq  13  were  the  more  valid  equa¬ 
tion,  the  slope  would  be  -3  f°r  Z  =  + 2.  Thus,  it  appears 
that  the  slope  of  -2.9  obtained  from  the  plot  in  Figure  27 
is  indicative  of  an  effective  charge  of  between  +2  and  +3 
on  the  P  species. 

The  heat  of  ionization  can  be  obtained  from  a  linear 
plot  of  log  Kr  vs.  l/T ;  the  slope  of  such  a  plot  is 
-4H°/2.303R.  The  ionic  strength  dependence  of  Kca  at  25.0° 
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Figure  26.  Plot  of  log  Kci  vs.  0.509(^ — ~  0.20ju) 
showing  the  ionic  strength  dependence  of  Kci  at  25*0°. 
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Figure  27*  Plot  of  log  Kc2  vs.  0.509(q — r~~  ~  0.20u) 

i  +  u 

showing  the  ionic  strength  dependence  of  Kc2  at  25»0°. 
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is  well  enough  defined  that  the  intercept  of  the  straight 
line  in  Figure  27  can  be  used  to  obtain  a  value  of  Kts  = 
(2.14*0.25 )  x  10” 11  M.  Values  of  KTa  at  18.0°  and  35.0°, 
calculated  from  eq  11,  using  a  value  of  -2.9  for  the  slope 
and  with  appropriate  correction  for  the  change  in  the 
Debye-Htickel  constant  with  temperature,  are  (1 .58*0. 10 ) x  10”1 1 
M  and  (3 • 33*0° 87)  x  10”11  M,  respectively.  The  logs  of 

'  r  ' 

these  values  vs.  1/T  are  plotted  in  Figure  28 j  from  linear 
least-squares  analysis,  AHa  is  found  to  equal  7° 8*0.1 
kcal/mole.  From  the  value  of  AH^j  additional  thermodynamic 
quantities  can  be  calculated  for  the  ionization  of  isozyme 
2  (type  C)  at  25«0°s 

AGa  =  -RTlnKTa  =  14.6*0.1  kcal/mole 

ASa  =  =  -22.6*1.0  e.u. 

These  values  may  be  compared  to  the  corresponding  thermo¬ 
dynamic  quantities  for  other  ionizations  of  water  in  differ¬ 
ent  environments,  shown  in  Table  X.  As  can  be  seen  in  that 
table,  the  thermodynamic  constants  for  the  ionization  of 
water  in  the  sixth  position  of  the  heme  proteins  are  not 

too  different  from  one  another,  but  as  a  group  they  differ 

/  \ 3  + 

from  either  the  ionization  of  water  alone  or  on  Fe(OHajs  ? 
showing  the  considerable  effect  of  the  protein-heme  struc¬ 
ture  on  the  nature  of  the  hydrogen- oxygen  bond  in  the  water 


molecule 


log  (Kt2  x  1011) 
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Figure  28.  Plot  of  log  Kt3  vs.  1/T  for  experiments  per- 

O  . 

formed  at  p.  =  0.11.  The  slope  is  -  H2/2.303R« 
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CHAPTER  ,6.  STUDIES  ON  HEME  PROTEINS  USING  THE  NMR  HALIDE 

ION  PROBE  TECHNIQUE 


The  use  of  halide  ions  as  probes  for  nuclear  magnetic 
resonance  studies  of  proteins  has  been  the  subject  of  sev- 
eral  recent  papers1 03"1 0 7 e  This  chapter  is  a  report  on  the 
use  of  the  halide  ion  probe  technique  to  study  four  heme 
proteins:  bovine  hemoglobin,  sperm  whale  myoglobin,  bovine 
liver  catalase,  and  horseradish  peroxidase* « 

Theory 

Nuclei  which  have  a  quadrupole  moment  can  interact 
with  a  fluctuating  electric  field  gradient  to  provide  an 
effective  nuclear  relaxation  mechanismo  For  3BC1,  which 
has  a  nucleus  of  spin  3/2  9  the  contribution  to  the  nuclear 
resonance  line  width  from  quadrupole  relaxation  is108: 

AV  =  ~(eaqQ)3  7^  (1) 

where  av  is  the  line  width  in  Hz  at  half-height;  the 

«»  UBD  CMTl  CB»  OD  CBS  Cl a  CZ)  CX3  Cii  GXJ  OB  OB  UB>  CKJ  CD  IB  LSD  OD  dO  IB  IQ  C3  CD  LSD  UD  OB  CSJ  OC1  OJ  UUI  CCJ  CX.»  CMJ  OBJ  CSI  (W>  CD  '  K3  US) 

"^Abbreviations  used  in  this  chapter:  Hb,  f  errihemoglobin ; 
Mb,  ferrimyoglobin ;  Cat,  catalase;  PCMB,  p-chloromercuri- 
benzoate;  P-S-Hg-  ,  protein  with  sulfhydryl  proton  replaced 
by  Hg3^  j  a^)obs  ,  total  observed  3  BC1  linewidth;  total 

3BC1  linewidth  minus  the  contribution  from  the  NaCl  present 
in  the  solution;  Hz,  Hertz,  or  cycles/sec 0 
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quantity  ( esqQ ) ,  known  as  the  quadrupole  coupling  constant, 
consists  of  e,  the  charge  on  the  electron,  q,  the  electric 
field  gradient  at  the  nucleus,  and  Q,  the  quadrupole  mo¬ 
ment  of  the  nucleus j  ^  is  the  rotational  correlation  time. 
A  large  range  of  linewidths  is  possible,  since  both  lc  and 
the  quadrupole  coupling  constant  can  vary  greatly  for  dif¬ 
ferent  molecules o  Dilute  aqueous  solutions  of  chloride  ion 
give  linewidths  of  about  15  Hz,  a  result  of  the  fact  that 
the  chloride  ion  is  essentially  symmetrically  solvated, 
which  produces  a  field  gradient  of  almost  zero  at  the  nu¬ 
cleus  o  If  the  chlorine  atom  is  involved  in  covalent  bond- 
ing,  a  large  value  of  the  coupling  constant,  and  hence,  a 
large  linewidth  results «  For  example,  the  linewidth  for 
CC14  would  be  of  the  order  of  14<>5  kHz1 03  0 

If  there  are  several  sites  available  for  the  chlorine, 
and  if  the  exchange  between  them  is  slow,  the  resulting 
spectrum  will  contain  lines  characteristic  of  the  different 
environments,  with  the  areas  proportional  to  the  concentra¬ 
tions  of  the  various  sites 0  If  the  exchange  between  the 
sites  is  fast  compared  to  the  reciprocal  of  the  width  of 
the  broadest  line,  a  single  composite  line  is  observed  with 
the  observed  linewidth,  ,  given  by? 

A  VMS  =  <5l  V>  P;  (  2  ) 

where  aV(-  is  the  natural  linewidth  of  3  6C1  in  the  i  th 
environment  and  P;  is  the  probability  that  the  chlorine  is 
at  the  i th  site0 
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Heme  proteins  used  in  this  study  were?  bovine  Hb 
(2X  ciystoj  A  grade*  lot  64385)  from  Calbiochem?  sperm 
whale  Mb  from  Mann  Research  Laboratories  $  analytical  reagent 
grade  HRP  (RZ  =  2c9®  lot  6485206)  and  purified  grade  HRP 
(RZ  a  O08)  from  Boehringer-Mannheim  Corpo ,  No  Y0 §  and  bovine 
liver  catalase  from  Pierce  Chemicals  (RZ  =  0d)  and  from 
Boehringer-Mannheim  Corpo*  N0  Yo  (analytical  reagent  grade* 
RZ  -  0u9s!  lot  6028247)o  Protein  solutions  were  dialyzed 
against  water  and  centrifuged  before  use.  All  other  chem¬ 
icals  used  in  this  study  were  reagent  grade  and  were  used 
without  further  purification 0  Concentrations  of  the  heme 
proteins  were  determined  from  absorbance  measurements  using 
either  a  Beckman  DU  or  a  Cary  Model  14  spectrophotometer, 
along  with  literature  values  for  the  molar  absorptivity 
coefficients  6  *  56 e  Measurements  of  pH  were  made  with  a 
Beckman  expanded  scale  pH  meter 0 

The  catalase  from  Pierce  Chemicals  was  further  purified 
by  ammonium  sulfate  precipitation  (50  °/o  saturation)  and 
dialysis,  followed  by  passage  through  a  2^5  cm  x  50  cm  col- 
umn  of  G-200  Sephadex  gel o  The  catalase  was  eluted  using 
pH  60 4  phosphate  buffer  of  ionic  strength  equal  to  0.1. 

The  purest  fraction  resulting  from  a  single  passage  through 
the  column  had  a  RZ  value  of  0o5« 

The  apoprotein  of  HRP  was  isolated  by  cleaving  analyt¬ 
ical  reagent  grade  HRP  at  pH  2 05  using  Yonetani's 
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modific.ation109  of  Teale8 s  2-butanone  method11 °0  Spectro- 
photometric  measurements  at  403  mp,  the  maximum  of  the 
Soret  peak  of  HRP,  indicated  that  less  than  2 <>5  °/0  of  the 
native  HRP  was  present  in  the  isolated  apo-HRP® 

PCMB  titrations  of  analytical  reagent  grade  catalase 
and  analytical  reagent  grade  HRP  were  carried  out  at  pH  6 »0 
using  the  method  of  Boyer111 0  Samples  were  allowed  to  stand 
at  room  temperature  for  24  hours  to  ensure  that  the  reaction 
was  complete  before  measurements  were  made  on  the  Cary 
Model  14 o 

The  3SC1  spectra  were  obtained  at  5°9  MHz  using  a 
Varian  VF16  wide  line  NMR  spectrometer  equipped  with  a 
V4012A  high  resolution  12n  magneto  A  Varian  V3521  inte- 

i 

grator  or  a  Princeton  Applied  Research  JB-4  lock-in  amp¬ 
lifier  were  used  for  field  modulation  and  detection 0  Line- 
widths  were  calibrated  from  sidebands  of  a  saturated  aqueous 
NaCl  solution o  A  Hewlett-Packard  24lA  oscillator  supplied 
the  sideband  frequency,  which  was  measured  using  a  Hewlett- 
Packard  5210  counter o  All  spectra  were  obtained  at  room 
temperature  using  5  ml  of  liquid  in  a  sample  tube  of  15  mm 
diameter®  Linewidths  reported  in  this  chapter  were  measured 
as  the  distance  between  the  turning  points  of  dispersion 
mode  spectra  and  are  the  average  of  at  least  five  spectra. 
Error  limits  reported  are  mean  deviations®  Representative 
spectra  are  shown  in  Figure  29  *> 

For  each  experiment,  spectra  were  taken  of  a  NaCl 


. 
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Figure  29.  Representative  NMR  spectra.  A.  Saturated  NaCl 
with  75  Hz  sidebands.  B.  3*^  x  10” 6  M  Hb  in  0.5  M  NaCl; 
linewidth  is  15  Hz.  C.  Broad  signal  with  a  linewidth  of 
65  Hz. 
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solution  of  a  molarity  equal  to  that  in  the  experimental 
solution o  In  the  case  of  the  experiments  involving  titra¬ 
tion  with  HgClSs  a  solution  of  protein  was  prepared  which 
was  either  0o5  M  or  loO  M  in  NaCl  and  was  titrated  with  pi 
volumes  of  a  solution  of  HgCla  which  was  about  100  times 
more  concentrated  than  the  protein «  This  procedure  elim¬ 
inated  the  necessity  for  a  correction  due  to  increase  in 
total  volume,  since  the  maximum  increase  in  volume  due  to 
added  titrant  was  less  than  5  °/o° 

Results 

For  solutions  used  in  the  present  study,  eq  2  can  be 
written  as  s 

AVaas  =  ^Hcr)P(cr>  +  (3) 

in  which  the  contribution  to  the  total  linewidth  from  the 
chloride  ion  in  solution  is  separated  from  that  from  chloride 
in  the  other  possible  environments 0  Because  the  chloride 
ion  concentration  was  at  least  104  times  that  of  the  protein, 
P(cr)  is  close  to  unity  and  eq  3  can  be  simplified  to* 

AV*  =  A\igs  ~  A}jcr)  =  (4) 

Further,  the  assumption  that  P(Cr)  ^  1  implies  that  the 

P'  s  are  independent  of  each  other 0  All  linewidths  obtained 
from  protein  solutions  will  henceforth  be  expressed  in  terms 
of  AV>\  i c e0 ,  the  linewidth  of  the  experimental  solution 
minus  the  linewidth  of  a  solution  containing  only  NaCl  of 
the  same  concentration  as  that  in  the  experimental  solution* 
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The  36C!  linewidth  of  a  solution  of  either  0*5  M  NaCl 
or  1  M  NaCl  at  pH  6*5  was  found  to  be  14  ±  2  Hz,  which  was 
not  changed  by  lowering  the  pH  to  2»0.  Addition  of  HgCla 
to  a  solution  of  NaCl  broadened  the  line  linearly?  a  0«5  M 
NaCl  solution  which  was  2  x  10“4  M  with  respect  to  HgCla 
gave  a  linewidth  of  25  ±  2  Hz* 

This  study  resulted  in  the  first  observation ,  by  the 
NMR  halide  ion  probe  technique,  of  direct  chloride  binding 
to  a  protein o  Figure,  30  shows  the  effect  of  pH  on  aV  c  for 
solutions  of  7  0  x  10“ 6  M  analytical  reagent  grade  HRP  and 
apo-HRP9  both  1  M  in  NaClo  HRP  precipitated  below  pH  3  and 
apo-HRP  below  pH  3°5»  thus  setting  the  lower  limits  of  pH 
at  which  the  study  could  be  carried  out.  Significant  broad¬ 
ening  due  to  chloride  ion  binding  to  catalase  was  observed 
for  solutions  of  catalase  obtained  from  Pierce  Chemicals* 

A  4d  x  10“ 6  M  solution  of  catalase  at  pH  6*5  which  was 
0*5  M  in  NaCl  gave  a  value  of  =50  *  5  The  amount 

of  broadening  was  the  same  for  solutions  of  catalase  from 
the  same  source,  but  which  had  been  further  purified*  The 
broadening  was  less  by  an  order  of  magnitude  for  a  solution 
of  the  analytical  reagent  grade  catalase?  for  2*8  x  10“ 6  M 
analytical  reagent  grade  catalase  at  pH  6*8  and  which  was 
0*5  M  in  NaCl,  =  6  ±  1  Hz* 

Curves  showing  the  titrations  of  analytical  reagent 
grade  HRP,  Hb,  and  analytical  reagent  grade  catalase  with 
HgCla,  a  reagent  which  is  quite  specific  for  sulfhydryl 


. 
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PH 


Figure  30*  Plot  of  line  broadening  vs.  pH  for  a  solu¬ 
tion  of  7.3  x  10“ 5  M  HRP  (o)  and  7.3  x  10“ 6  M  apo-HRP 
(x),  both  1  M  in  NaCl . 
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groups1 2,3  *  are  shown  in  Figures  31,32,  and  33,  with  the  re¬ 
sults  plotted  as  &))  J  vss  [HgCla ]/[protein] „  The  titration 
curves  of  HRP  and  catalase  indicate  the  presence  of  one  and 
two  sulfhydryl  groups 9  respectively j  the  titration  curve  of 
Hb  is  not  as  well  defined,  appearing  to  break  between 
[HgCla ]/£Hb]]  ratios  of  one  and  two0  The  HgCla  titration 
curve  for  purified  grade  HRP  was  identical  to  that  of  ana¬ 
lytical  reagent  grade  HRPs  Spectrophotometric  titration  of 
analytical  reagent  grade  HRP  and  analytical  reagent  grade 
catalase  with  PC MB  indicated  the  presence  of  l«0±0o2  and 
4o5±0c5  sulfhydryl  groups,  respectively 0 

The  results  of  the  titration  of  4o 46  x  10"s  M  Mb  at  pH 
7*9  an(3  [Cl°°]  =  1M  with  HgCla  are  shown  in  Figure  3^°  The 
dashed  line  in  that  figure  shows  the  broadening  expected 
from  the  addition  of  HgCla  alone,  i«e°  ,  if  there  were  no 
interaction  between  the  HgCla  and  the  myoglobin 0  For  a 
given  amount  of  HgCla  added  to  the  Mb  solution,  the  line  — 
width  was  found  to  be  pH  dependents  This  is  illustrated 
in  Figure  35?  which  is  a  plot  of  A))'1  vs o  pH  for  a  sample 
containing  1  M  NaCl,  4o5  x  10“ 6  M  Mb,  and  hi  x  10“° 4  M 
HgCla  o  The  broadening  expected  if  there  were  no  interac¬ 
tion  between  Mb  and  HgCla  is  again  shown  by  a  dashed  lines 
There  appears  to  be  no  break  indicating  the  presence  of  an 
end  point  in  Figure  3^»  another  titration  carried  out  to  a 
ratio  of  [HgCla]/[Mb]  of  greater  than  10  also  indicated  no 
breaks  The  reversibility  of  the  reaction  between  HgCla 
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Figure  31.  Curve  showing  the  titration  of  5*8 
M  HRP  with  HgCl2;  pH  =  7.0,  [Cl~]  =  0.5  M. 


-125- 


[HgCla] 
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Figure  J2 •  Curve  showing  the  titration  of  3»^  x  10~5 
M  Hb  with  HgCl3  ;  pH  =  6. 3»  [Cl“]  =  0.5  M. 
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Figure  33.  Curve  showing  the  titration  of  1.0  x  10  5 
M  catalase  with  HgCl2 ;  pH  =  6.6,  [Cl~]  =  1  M. 
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Figure  3^«  Curve  showing  the  titration  of  4.46  x  10~5 
M  Mb  with  HgCl2;  pH  =  7»9j  [Cl~]  =  1  M.  The  dashed 
line  shows  the  broadening  expected  if  there  were  no 
interaction  between  the  Mb  and  the  HgCl2. 
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Figure  35*  Plot  of  line  broadening  vs.  pH  for  a  solu¬ 
tion  of  4.5  x  1CT5  M  Mb;  [HgCls]  =  1.1  x  1CT4  M; 

[Cl“]  =  1  M.  The  dashed  line  shows  the  broadening 
expected  if  there  were  no  interaction  between  the  Mb 
and  the  HgCl2. 
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and  Mb  was  demonstrated  by  lowering  the  pH  of  a  solution 
containing  1  M  NaCl,  2d  x  10" 6  M  Mb,  and  2»1  x  10” 4  M 
HgGla  from  pH  9° 2  to  5»8«.  This  resulted  in  1  decreasing 
from  56  ±  5  Hz  to  11  ±1  Hz9  the  value  expected  from  the 
presence  of  unreacted  HgClas 

The  reactivity  of  HgCla  with  methionine  and  imidazole 
was  investigated  at  pH  values  of  7  and  10  by  observing  the 
linewidths  of  solutions  of  1  M  NaCl  containing  3®^  x  10” 4  M 
HgCls  with  and  without  the  presence  of  either  4  x  10“ 4  M 
methionine  or  4  x  10“ 4  M  imidazole 0  No  evidence  of  inter¬ 
action  between  HgCla  and  either  compound  was  observed 0 

Discussion 

The  observed  binding  of  chloride  ions  to  catalase  and 
analytical  reagent  grade  HRP  demonstrates  the  potential 
usefulness  of  the  halide  ion  probe  technique  as  a  tool  for 
observing  protein-chloride  interaction 0  Since  it  has  been 
reported  that  the  Soret  absorption  peak  of  HRP  undergoes  a 
spectral  change  at  low  pH  in  the  presence  of  chloride,  a 
result  which  was  attributed  to  the  binding  of  chloride  to 
the  heme  ferric  iron  of  HRP3  0 ,  one  might  expect  that  the 
observed  line  broadening  was  due  to  chloride  exchange  at 
the  sixth  coordination  position  of  the  iron,  and  that  the 
observed  pH  dependence  is  correlated  with  the  pH  dependence 
of  the  equilibrium  constant  for  HRP-chioride  formation 0 
However,  this  explanation  cannot  be  correct,  since  the  pH 
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dependence  of  chloride  binding  to  apo-HRP  is  similar  to 
that  of  native  HRPS  as  can  be  seen  in  Figure  30°  The 
results  of  this  study  do  not  rule  out  the  possibility  that 
chloride  binds  to  the  heme2  since  for  the  heme-chloride 
site  could  simply  be  too  small  to  result  in  a  contribution 
to  the  total  linewidth8  as  expressed  in  eq  2°  The  pH  de¬ 
pendence  of  the  chloride  binding  could  be  explained  if  the 
chloride  were  binding  electrostatically  and  nonspecif ically 
in  regions  of  appreciable  positive  charge^  since  the  net 
positive  charge  on  the  protein  would  increase  as  the  pH  was 
lowered  due  to  protonation  of  various  groups °  This  hypoth¬ 
esis  is  substantiated  by  the  fact  that  the  amount  of  broad¬ 
ening  at  a  given  pH  is  greater  for  apo-HRP  than  it  is  for 
native  HRP0  An  ORD  study  has  shown  that  as  the  heme  is 
split  from  the  protein  of  HRPS  the  protein  becomes  less 
helical  (Chapter  4  of  this  thesis )*  thus  the  greater 
chloride  binding  for  apo=HRP  could  be  a  result  of  the  par¬ 
tial  unfolding  of  the  protein  and  the  exposure  of  more 
binding  sites °  One  can  also  see  in  Figure  30  that  the 
amount  of  line  broadening  for  HRP  increased  apprecially 
below  pH  3° 4S  which  might  be  due  to  the  rapid  splitting  of 
the  heme  from  the  protein  of  HRP300 

The  results  of  the  study  of  chloride  binding  to  cata¬ 
lase  are  difficult  to  explain ,  since  they  varied  with  the 
source  of  the  catalase °  There  are  two  possibilities  which 
could  account  for  the  observed  results °  The  first  is  that 
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there  is  some  impurity  capable  of  binding  chloride  present 
in  the  catalase  preparation  obtained  from  Pierce  Chemicals® 
If  this  is  the  case,  the  impurity  must  be  protein  material, 
since  it  was  not  removed  by  exhaustive  dialysis 0  However , 
further  purification  of  that  preparation  increased  the 
ratio  of  absorbances  of  the  maxima  of  the  heme  and  protein 
peaks  (RZ)  by  a  factor  of  five,  indicating  that  a  consider” 
able  amount  of  protein  was  removed  during  the  further  puri- 
ficationo  Thus,  if  the  observed  results  are  due  to  a  pro¬ 
tein  impurity,  the  protein  must  have  a  molecular  weight 
similar  to  that  of  catalase  since  the  same  line  broadening 
was  observed  before  and  after  additional  purification „  An 
alternative  explanation  would  be  that  the  isolation  proce¬ 
dure  used  by  Pierce  Chemicals  resulted  in  a  change  in  the 
protein  such  that  groups  capable  of  binding  chloride  were 
exposed a 

This  investigation  also  included  studies  of  the  behav- 
ior  of  Hga+  as  a  sulfhydryl  reagent  when  the  halide  ion 
probe  method  is  used  as  a  means  of  detection »  The  titra¬ 
tion  of  HRP  with  HgCls  has  a  well-defined  end  point  indicat¬ 
ing  one  sulfhydryl  group,  as  is  shown  in  Figure  31°  The 
positive  slope  beyond  the  end  point  in  Figure  31  is  due  to 
the  contribution  from  uncomplexed  HgClao  The  presence  of 
one  sulfhydryl  group  in  HRP  was  confirmed  by  the  spectro- 
photometric  titration  using  PCMBo  No  previous  quantitative 
study  has  been  made  on  the  sulfhydryl  groups  of  HRP,  but 
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qualitative  amino  acid  analysis  has  shown  the  presence  of 
cysteine1 6 o  Although  purified  grade  HRP  contains  a  greater 
amount  of  protein  material  than  does  the  analytical  reagent 
grade  HRP9  as  evidenced  by  its  considerably  lower  RZ  value, 
results  from  the  two  protein  preparations  were  the  same, 
indicating  the  presence  of  no  additional  reactive  sulfhydryl 
groups  in  the  protein  impurities  present » 

It  had  previously  been  reported  that  good  results  were 
obtained  using  this  method  for  the  titration  of  equine  and 
human  hemoglobin1 04 0  However,  the  HgCla  titration  curve 
for  bovine  Hb,  shown  in  Figure  32,  is  difficult  to  inter¬ 
pret  o  It  is  obvious  that  there  is  at  least  one  sulfhydryl 
group,  and  probable  that  there  is  a  second,  since  the  titra¬ 
tion  curve  has  a  slope  characteristic  of  HgCla  alone  only 
after  a  [HgCla ]/[Hb]  ratio  of  two®  Direct  comparison  of 
this  curve  with  those  obtained  by  Stengle  and 
Baldeschwieler 1 04  cannot  be  made  for  two  reasons o  First, 
their  data  is  for  equine  and  human  Hb,  while  this  study  was 
carried  out  using  bovine  Hb,  and  although  they  found  the 
same  number  of  sulfhydryl  groups  for  both  of  those  hemo¬ 
globins  (two),  there  could  be  a  species  difference  between 
those  two  types  and  bovine  Hb®  Second,  their  published 
curves  show  points  only  at  integral  [HgCla ]/[Hb]  ratios, 
so  any  anomalous  behavior  might  not  show  up  due  to  a  lack 
of  sufficient  points®  The  simplest  explanation  for  the 
behavior  shown  in  Figure  32  is  that  there  are  two  reactive 
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sulfhydryl  groups  in  bovine  Hb,  but  that  they  are  in 
environments  of  different  accessibility*  Thus,  if  the 
second  site  is  less  accessible  and  has  a  smaller  7c  than 
the  first,  the  slope  of  the  titration  curve  of  the  second 
site  would  be  less  than  that  of  the  first,  as  predicted  by 
eq  1. 

The  results  of  the  HgCla  titration  of  catalase,  which 
indicate  that  there  are  two  reactive  sulfhydryl  groups, 
disagree  with  the  results  of  the  PCMB  titration  which 

indicates  4*5::fc0*5  reactive  sulfhydryl  groups,  and  both  of 
these  disagree  with  the  findings  of  Samejima  and  Yang113, 
who  found  six  reactive  groups  in  catalase  by  titration  with 
PCMBo  The  disagreement  between  the  results  obtained  from 
the  two  methods  could  be  due  to  the  nature  of  the  action  of 
the  two  different  sulfhydryl  reagents*  HgCla  reacts  rapidly 
and  perhaps  only  with  the  most  accessible  sulfhydryl  groups* 
PCMB  reacts  slowly  with  catalase,  and  it  has  been  suggested 
that  in  such  cases,  some  type  of  bond  has  to  be  broken 
before  the  sulfhydryl  groups  can  react113,  thus  the  PCMB 
could  be  reacting  with  additional,  less  reactive  sulfhydryl 
groups*  The  slow  reaction  of  PCMB  with  catalase  could  also 
explain  the  difference  between  the  results  of  this  study 
and  those  obtained  by  Samejima  and  Yang*  Those  authors 
did  not  tell  under  what  conditions  they  found  the  six 
titratable  groups,  and  since  they  found  an  additional 
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ten  sulfhydryl  groups  upon  denaturation  of  the  catalase, 
it  is  possible  that  there  exists  in  catalase  sulfhydryl 
groups  of  many  degrees  of  reactivity,  and  that  Samejima  and 
YangTs  conditions  for  the  PCMB  titration  were  such  that  more 
groups  could  react  than  under  the  conditions  used  in  this 
study.  However,  as  pointed  out  in  the  review  by  Cecil  and 
McPheella,  results  of  determinations  of  the  number  of 
sulfhydryl  groups  of  proteins  often  differ  depending  upon 
the  sulfhydryl  reagent  used,  and  also  that  different  inves¬ 
tigators  frequently  obtain  conflicting  results  using  the 
same  sulfhydryl  reagent. 

The  results  of  the  titration  of  Mb  with  HgCla  are  quite 
surprising  because  Mb  is  known  to  contain  no  cysteine  or 
cystine114.  Since  mercury  is  a  member  of  the  same  group 
as  zinc,  and  since  Zna  +  is  known  to  bind  to  imidazole  res¬ 
idues  in  Mb116'116,  there  is  a  possibility  that  Hga  +  might 
also  react  with  imidazole.  Likewise,  since  methionine  con¬ 
tains  sulfur,  which  although  incapable  of  forming  a  complex 
of  the  type  R-S-Hg-Cl,  might  nevertheless  be  capable  of 
forming  a  complex  with  mercury  through  one  of  its  unshared 
electron  pairs.  However,  some  evidence  against  either 
imidazole  or  methionine  residues  providing  binding  sites  is 
furnished  by  the  lack  of  reactivity  of  HgCla  with  either 
compound.  The  lack  of  any  break  in  Figure  3^  indicates 
that  either  (a)  the  number  of  binding  sites  is  large,  or 
(b)  the  equilibrium  constant  for  the  formation  of  the 
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mercury-protein  complex  is  very  small. 

It  is  interesting  to  note  that  Figure  35  looks  much 
like  a  mirror  image  of  Figure  30.  Therefore,  one  might 
postulate  that  the  observed  pH  dependence  shown  in  Figure 
35  is  a  result  of  Hgs+  binding  to  anionic  sites  on  the  pro¬ 
tein,  and  that  the  number  of  anionic  sites  increases  as  the 
pH  is  raised. 

This  study  has  shown  the  halide  ion  probe  technique 
to  be  a  very  sensitive  method  for  detecting  the  binding  of 
Cl"  or  Hga+  to  proteins.  When  used  as  a  means  of  detection 
for  titrations  with  HgCla  to  determine  the  sulfhydryl  con¬ 
tent  of  proteins,  this  method  has  the  advantages  of  speed 
and  simplicity  and  offers  the  possibility  of  distinguishing 
sulfhydryl  sites  of  different  accessibility.  Judgment  of 
the  relative  merit  of  this  means  of  determining  the  sulf¬ 
hydryl  content  of  proteins,  as  compared  with  older  tech¬ 
niques,  must  await  further  studies. 
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APPENDIX  I,  COMPUTER  PROGRAMS 


This  appendix  contains  user  information  and  listings 
of  the  two  basic  computer  programs  used  throughout  the 
course  of  the  investigations  reported  in  this  thesis « 

Nonlinear  Least-Squares 

A  nonlinear  least-squares  program  is  invaluable  for  the 
fitting  of  data  to  equations,  such  as  eq  35  of  Chapter  3» 
which  contain  a  large  number  of  adjustable  parameters.  A 
nonlinear  least-squares  program  was  obtained  from  the  IBM 
Share  Library  ( SDA  3094)  and  was  modified  to  use  double 
precision  throughout  and  to  run  on  the  University  of  Alberta 
IBM  36O  computer. 

The  program  contains  three  subroutines  for  which  the 
user  must  supply  appropriate  Fortran  statements.  Immedi¬ 
ately  after  the  data  is  read  into  the  computer,  the  program 
calls  subroutine  SUBZ,  which  may  either  be  left  blank  or 
used  at  the  option  of  the  programer.  This  subroutine  can 
be  used  to  calculate  constants  needed  in  the  model  equation, 
or  to  perform  auxiliary  calculations,  such  as  conversion  of 
pH  to  [H+],  before  the  main  program  is  executed. 

Subroutine  F  CODE  contains  the  model  equation  to  which 
the  data  is  to  be  fit®  The  model  must  be  of  the  form 
F  =  f (x( 1 , 1 ) ,  X(l,2) , . » 0  j  B ( 1 ) ,  B(2),00o),  where  F  predicts 
the  value  of  the  dependent  variable  Y(l);  X(l,l),  X(l,2),.o» 


' 
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are  the  independent  variables  §  and  B(l),  B(2),0..  are  the 
adjustable  parameters  or  the  unknown  constants  of  the  model 
equation. 

The  subroutine  P  CODE  evaluates  the  analytic  deriva¬ 
tives  of  the  model  function  for  the  same  point  just  pro¬ 
cessed  by  F  coding.  If  analytic  derivatives  are  not  sup¬ 
plied,  i.eo ,  if  P  CODE  is  left  blank  (and  ISW2  is  set  equal 
to  1  in  input  item  number  2),  the  program  will  estimate 
values  of  BF;,/^>B(l)  using  finite-difference  approximations. 
Initial  guesses  for  the  parameters  must  be  different  from 
zero  if  estimated  derivatives  are  used. 

Once  the  subroutines  have  been  set  up,  the  data  cards 
are  prepared  as  follows  s 

Input  Item  Fortran 

Number  Label  Format  _ Comment _ 


1 

N 

13 

no.  of  data  points  (Nmr*=500 

K 

13 

total  no.  of  parameters 
(KMax=  25) 

IP 

13 

must  be  zero 

M 

13 

no.  of  independent  variables 
(Mh«,=  10) 

2 

ISW1 

13 

must  be  zero 

ISW2 

13 

=.0 ,  analytic  derivatives 
=1 ,  estimated  derivatives 

isw3 

13 

=0,  abbreviated  print  out 

=1 ,  detailed  print  out 

Detailed  print  out  gives  a  complete  statement  of  observed 
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and  predicted  values  of  the  dependent  variables  for  each 
new  choice  of  trial  values  for  the  parameters 0 


Input  Item 
Number 

Fortran 

Label 

Format 

Comment 

3 

B(l) 

DIO.  0 

initial  guesses  for  param' 
eters,  7  per  card 

B(2) 

DIO .  0 

etc. 

DIO  0  0 

4.  This  item  consists  of  N  sub-items <>  Each  sub-item  is 
the  input  data  for  one  observations,  punched  accord¬ 
ing  to  the  format  8D10.0.  The  sequence  of  variables 
in  a  sub- item  must  be  Y(l),  X(I,l),  X(  1, 2  ),«..> 
through  X(I,M).  Each  Y(l)  must  begin  on  a  new  card* 

5c  Any  case  data  read  in  from  subroutine  SUBZ  should 
go  here. 

60  Sequential  cases  may  be  stacked  by  repeating  items  1 
through  5°  A  blank  card  after  the  last  case  will 
cause  a  normal  stop  without  an  error  message 0 

Weighting  Modification  of  the  Nonlinear  Least-Squares 

Program 

If  the  observed  values  of  the  dependent  variable  vary 
by  more  than  one  order  of  magnitudes  a  ” skewed w  fit  to  the 
data  may  result  if  the  program  is  used  in  its  original  form. 
This  problem  can  best  be  understood  by  looking  at  an  exam¬ 
ple.  Suppose  a  series  of  apparent  rate  constants  exist 
which  are  to  be  fit  by  a  model  equation  and  which  vary  from 
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k(l)  =  loO  x  104  to  k(lO)  =  loO  x  10s,  with  the  experimental 
errors  of  all  the  points  equal  at  about  =fcl0°/o°  If  the 
model  assumes  trial  values  of  the  parameters  such  that  the 
predicted  rate  constants  are  k(l)  =  Id  x  106  and  k(10)  = 

Id  x  106,  the  difference  between  the  observed  and  predic¬ 
ted  values  in  each  case  is  L0  x  106,  It  is  the  sum  of  the 
squares  of  these  residuals  which  the  program  attempts  to 
minimize 0  However,  in  thq  case  above,  it  can  be  seen  that 
the  predicted  value  of  k(l)  is  10  times  the  observed  value, 
but  the  predicted  value  of  k(10)  agrees  within  10  °/o  of 
the  observed  value 0  A  simple  remedy  for  this  inequitable 
weighting  is  the  following  DO  loop  which  is  to  be  inserted 
into  subroutine  SUBZ0  For  this  example,  it  is  assumed  that 
the  Nth  observed  dependent  variable  is  the  largest  and  that 
there  is  only  one  independent  variable.. 

DO  1  1=1, N 
X(I,2)  =  y(n)/y(i) 

1  Y(l)  =  Y(l)*X(l,2) 

By  this  means,  we  have  created  a  second  independent  vari¬ 
able  which  is  in  reality  a  weighting  factors 

Now,  immediately  following  the  statement  defining  the 
model  in  subroutine  F  CODE  should  go  the  statements 
F  =  F*X(I,2) 

Thus,  both  the  observed  and  the  predicted  values  have  been 
multiplied  by  the  weighting  factor  (or  in  this  case,  more 
properly,  a  normalizing  factor )»  This  weighting  method  is 
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also  applicable  if  the  user  wishes  to  weight  the  data 
points  according  to  their  experimental  uncertainties,  etc. 
One  then  merely  needs  to  have  the  weighting  factor  read  in 
as  an  independent  variable  and  to  use  the  statements  de- 
scribed  above,  but  omitting  the  middle  card  of  the  DO  loop. 

A  listing  of  the  nonlinear  least-squares  program  and 
its  subroutines  is  found  on  the  following  pages. 

Polynomial  Least- Squares 

A  polynomial  least-squares  program,  obtained  from 
Dr.  J.  Ac  Plambeck,  was  used  for  all  linear  least-squares 
analyses  which  were  performed.  This  program  is  found  fol¬ 
lowing  the  listing  of  the  nonlinear  least-squares  program; 
the  listing  contains  user  instructions  in  comment  state¬ 
ments  at  the  beginning. 
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APPENDIX  II 

PURIFICATION  OF  HRP  USING  SEPHADEX  GELS 


The  purity  of  HRP  is  commonly  expressed  in  terms  of  a 
Reinheitszahl ,  RZ,  or  purity  number,  PN0  This  quantity  is 
determined  by  taking  the  ratio  of  the  absorbance  of  the 
maximum  of  the  Soret  peak  (403  mp)  to  that  of  the  maximum 
of  the  protein  peak  (280  mp)®  The  maximum  RZ  for  HRP  is 
3o0414a  HRP  is  commercially  available  from  Boehringer- 
Mannheim  Corp0 ,  No  Y0,  in  two  forms;  purified  grade,  at 
about  $45 *  00/g ,  with  an  RZ  of  about  O08;  and  analytical 
reagent  grade,  at  about  4l000o00/g,  with  an  RZ  ranging  from 
2*9  to  3°0o  This  study  was  undertaken  to  provide  HRP  of  an 
intermediate  purity  for  the  studies  described  in  Chapter  5 
and  to  determine  if  the  crude  HRP  could  be  economically 
and  easily  purified^  An  earlier  qualitative  study  by 
another  worker  in  this  laboratory,  Raymond  Segal,  indicated 
that  Sephadex  gels  would  provide  a  useful  medium  for  the 
purification  of  HRP 0 


Preparation  of  Gels  and  Columns 


Two  types  of  Sephadex  gels  were  used  in  the  purifica¬ 
tion  tests;  G-75s  a  molecular-sieve  with  a  molecular  weight 
fractionation  range  of  1000-50*000?  and  C-50  CK-Sephadex, 
which  is  a  carboxymethyl  derivative  of  Sephadex  that  func¬ 
tions  as  a  cation  exchanger 0  The  two  types  of  Sephadex 


' 


-174- 


were  prepared  for  use  in  the  following  manner:  3*5  9  of 

G-75  and  1.5  g  of  CM  were  each  allowed  to  swell  for  24 

1 

hours  in  about  15 0  ml  of  phospate  buffer  of  pH  6.2  and  p 
equal  to  0.1,  with  three  changes  of  supernatant  liquid  to 
ensure  complete  equilibration.  Approximately  50  ml  of  each 
type  of  swollen  gel  were  obtained. 

Chromatographic  columns  were  constructed  using  50  ml 
burets.  A  small  plug  of  glass  wool  was  placed  at  the  bot¬ 
tom  of  the  column  and  some  3  mm  diameter  glass  beads  were 
placed  on  top  of  the  glass  wool.  The  columns  were  filled 
by  pouring  a  slurry  of  the  gel  into  the  column.  After  a 
layer  of  gel  a  few  centimeters  thick  had  formed  at  the  bot¬ 
tom,  the  outlet  was  opened  to  aid  settling.  Each  column  was 
filled  until  it  had  a  bed  length  of  about  50  cm.  In  the 
case  of  the  CM  column,  the  bed  length  was  not  constant, 
since  the  size  of  the  gel  network  is  a  function  of  ionic 
strength.  The  columns  were  flushed  with  25  ml  of  buffer 
prior  to  addition  of  the  sample.  Because  the  gel  bed  com¬ 
pacted,  the  flow  through  the  column  was  quite  slow  even 
with  the  stopcocks  completely  open,  and  further,  it  was  not 
completely  constant  throughout  the  duration  of  the  separa¬ 
tion.  These  studies  were  conducted  with  the  stopcocks  of 
the  columns  completely  open,  with  flow  rates  of  the  order 
of  0.5  ml/min. 


- 
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Purification  and  Performance  of  the  Columns 

The  concentration  of  the  purified  grade  (RZ  =  O.77) 

HRP  was  determined  spectrophotometrically  to  be  7.6  x  10“4 
M  using  light  of  497  nip  and  a  molar  absorptivity  value8  of 
1.00  x  104  M”1  cm*"1.  One  ml  of  the  purified  grade  HRP  was 
added  to  each  column  and  when  the  solution  drained  to  the 
top  of  the  bed,  two  ml  of  pH  6.2  phosphate  buffer,  p  =  0.1, 
were  added.  The  buffer  did  not  acquire  any  brownish  color 
characteristic  of  the  HRP  solution,  which  indicated  that 
all  of  the  HRP  had  been  adsorbed  on  the  columns. 

The  HRP  was  eluted  through  the  G-75  column  entirely 
with  buffer  of  ionic  strength  0.1,  but  with  the  CM  column, 
the  ionic  strength  was  gradually  increased  from  0.1  to  0.2 
to  aid  the  separation.  Fractions  of  effluent  were  collected 
from  the  columns  and  the  concentration  of  HRP  in  each  frac¬ 
tion  was  determined  spectrophotometrically.  The  results  of 
the  single  column  purifications  are  shown  in  Tables  XI  and 
XII.  After  the  fourth  fraction  of  HRP  was  collected  from 
the  G-75  column,  the  column  showed  no  brown  coloration. 
However,  after  the  fifth  fraction  was  collected  from  the 
CM  column,  there  remained  a  small  brown  band  at  the  top  of 
the  column,  and  the  entire  column  had  a  slightly  brownish 
tinge.  The  top  band  started  to  move  only  after  the  ionic 
strength  of  the  buffer  was  increased  to  approximately  0.5 
by  addition  of  potassium  nitrate.  By  using  buffer  of  ionic 
strength  up  to  0.8,  the  CM  column  was  completely  cleaned  of 
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Table  XI 


Purification  of 

Crude  HRP  on  a 

g-75 

Sephadex  Column 

Fraction 

Ml 

Moles  x  108 

RZ 

°/o  Yield* 

1 

2.6 

5.0 

0.4 

6 . 6 

2 

5-0 

52.5 

1.6 

69.2 

3 

4.0 

14.8 

0.9 

19.5 

4 

3.4 

2.5 

0.2 

30 

98.6 

Table  XII 

Purification  of  Crude  HRP  on  a  CM-Sephadex  Column 


Fraction 

Ml 

Moles  x 

108 

RZ 

°/o  Yield* 

1 

2.7 

5.3 

1.9 

7.0 

2 

5 « 6 

34.7 

2.8 

45.8 

3 

2.9 

14.1 

2.7 

18.6 

4 

3.0 

4.4 

1.8 

5 « 8 

5 

4.6 

1.5 

0 . 6 

2.0 

6 

6 . 0 

3.0 

0.2 

3.9 

83.I 

*7.58  x  10" 

7  moles 

of  HRP, 

RZ 

=  0.77, 

were  added 

initially. 
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all  traces  of  brownish  material. 

After  dialysis  to  get  rid  of  KNOa  and  buffer,  the 
spectrum  of  fraction  number  six  (the  slow  fraction)  from 
the  CM  column  was  examined.  The  Soret  peak  of  this  fraction 
had  a  maximum  at  the  same  wavelength  (403  mp)  as  the  other 
fractions,  but  exhibited  a  more  gradual  decrease  of  absor¬ 
bance  with  decreased  wavelength.  The  higher  absorbance  on 
the  UV  side  of  the  peak  may  be  due  to  the  very  intense 
protein  peak  in  this  fraction  (RZ  =  0.2). 

In  order  to  ascertain  if  this  technique  can  be  used 
to  obtain  fractions  purer  than  those  shown  in  Tables  XI  and 
XII,  part  of  fraction  number  two  from  the  G»75  column  was 
passed  through  the  CM  column.  It  was  felt  that  the  two 
types  of  columns  used  in  series  might  be  more  effective 
than  a  second  passage  through  the  same  column.  1.16  x  10” 7 
moles  of  G-75  fraction  number  two  were  added  to  the  CM  col¬ 
umn  and  eluted  as  described  above;  the  results  are  shown 
in  Table  XIII. 

Table  XIII 

Further  Purification  of  G-75  Purified  HRP  on  a  CM  Column 


Fraction 

L  Ml 

Moles  x  108  RZ 

°/o  Yield* 

1 

2.6 

1.92  2.74 

16.5 

2 

2.6 

3.22  2.98 

27.8 

3 

2.6 

2.15  2.86 

18.5 

4 

2.6 

0.75  2.26 

6.5 

69.3 

*1 .16  x  10"7 

moles  of 

HRP,  RZ  =  1 .6,  were 

added  initially. 
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As  before,  a  small  band  was  left  at  the  top  of  the  column. 
This  band  was  not  collected  because  of  its  small  size.  If 
one  calculates  the  yield  from  the  original  crude  HRP,  frac¬ 
tion  number  two  of  the  doubly  purified  HRP  amounted  to 
19 • 2  °/o  of  the  original,  and  fraction  number  three  12.8  °/o» 
These  two  fractions  combined  represented  J2  °/o  yield  from 
the  original  purified  grade  HRP  used,  with  a  corresponding 
purification  from  RZ  =  0.7?  to  RZ  =  2.93* 

Although  no  attempt  was  made  to  maximize  either  purifi¬ 
cation  or  yield  through  variation  of  flow  rates,  column 
lengths,  etc.,  it  has  been  clearly  demonstrated  in  this 
study  that  HRP  can  be  purified  easily  with  good  yields 
using  Sephadex  gels. 
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